Theory of resonant multiphonon Raman scattering in graphene 
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We present a detailed calculation of intensities of two-phonon and four-phonon Raman peaks 
in graphene. Writing the low-energy hamiltonian of the interaction of electrons with the crystal 
vibrations and the electromagnetic field from pure symmetry considerations, we describe the system 
in terms of just a few independent coupling constants, considered to be parameters of the theory. 
The electron scattering rate is introduced phenomenologically as another parameter. 

The results of the calculation are used to extract information about these parameters from the 
experimentally measured Raman peak intensities. In particular, the Raman intensities are sensitive 
to the electron scattering rate, which is not easy to measure by other techniques. Also, the Raman 
intensities depend on electron-phonon coupling constants; to reproduce the experimental results, one 
has to take into account renormalization of these coupling constants by electron-electron interaction. 
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I. INTRODUCTION 

In the past decades, Raman spectroscopy^ techniques 
were successfully applied to carbon compounds, such as 
graphite (see Ref. 2 and references therein) and carbon 
nanotubes.^'^ Upon the discovery of graphene,^ Raman 
spectroscopy has proven to be a powerful tool to iden- 
tify the number of layers, structure, doping, disorder, 
and to characterize the phonons and electron-phonon 
coupling. ^"^^ So far, most of the attention was focused 
on the position and width of the Raman peaks. 

Here we present a detailed calculation of the intensities 
of the multiphonon Raman peaks in graphene. Raman 



scattering involves an electron-hole pair as an intermedi- 
ate state; we show that the multiphonon Raman peaks 
are strongly sensitive to the dynamics of this electron- 
hole pair. Thus, Raman scattering can be used as a tool 
to probe this dynamics. Writing the low-energy hamilto- 
nian of the interaction of electrons with the crystal vibra- 
tions and the electromagnetic field from pure symmetry 
considerations, we describe the system in terms of just a 
few independent coupling constants, considered to be pa- 
rameters of the theory. The electron scattering rate is in- 
troduced phenomenologically as another parameter. The 
results of the present calculation are used to extract in- 
formation about these parameters from the Raman peak 
intensities, measured experimentally. 

As shown below, the Raman intensities strongly de- 
pend on the electron scattering rate; moreover, the 
electron-phonon and electron-electron contributions to 
this rate can be separated. This is especially important 
as there are very few techniques giving experimental ac- 
cess to electron scattering rates, which, in turn, deter- 
mine the transport properties of graphene samples. Be- 
sides, the quasiclassical character of the process imposes 
a severe restriction on the electron and hole trajectories 
which can contribute to the two-phonon Raman scatter- 
ing: upon the phonon emission the electron and the hole 



J 



must be scattered backwards. This restriction results in a 
significant polarization memory: it is almost three times 
more probable for the scattered photon to have the same 
polarization as the incidend photon than to have the or- 
thogonal polarization. 

Also, the Raman intensities depend on electron- 
phonon coupling constants; to reproduce the experimen- 
tal results, one has to take into account renormalization 
of these coupling constants by electron-electron interac- 
tion. This renormalization is missed by local or semi- 
local approximations to the density-functional theory, 
typically used for the ab initio calculation of the coupling 
constants. 



A. Fully resonant processes 

Since graphene is a non-polar crystal, Raman scatter- 
ing involves electronic excitations as intermediate states: 
the electromagnetic field of the incident laser beam inter- 
acts primarily with the electronic subsystem, and emis- 
sion of phonons occurs due to electron-phonon interac- 
tion. The matrix element of the process can be schemat- 
ically represented as 
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Here \i) is the initial state of the process (the incident 
photon with a given frequency and polarization, and no 
excitations in the crystal), |/) is the final state (the emit- 
ted photon and n phonons left in the crystal), while Sk, 
fc = 0, . . . , n, label the intermediate states where no pho- 
tons are present, but an electron- hole pair is created in 
the crystal and k phonons have been emitted. Ei and Ek, 
k = 0, . . . , n are the energies of these states, and 27 is 
the inverse lifetime of the electron (hole) due to colli- 
sions. He-em and He-ph staud for the terms in the sys- 
tem hamiltonian describing interaction of electrons with 
the electromagnetic field and with phonons, respectively. 

In the calculations we do not include the phonon 
broadening, assuming the phonon states to have zero 
width. First, this approximation is consistent with the 
available experimental information: the phonon width is 
about 10 — 20 cm~^ « 2 — 3 meV at most,®'^° while the 
electronic broadening is at least an order of magnitude 
higher (see the discussion below. Sec. II C). Second, this 
approximation is irrelevant provided that we calculate 
the integrated intensities of the Raman peaks, since they 
are determined by the total spectral weight of the phonon 
state which does not depend on the phonon broadening. 



The photon wave vector is negligible, so momentum 
conservation requires that the sum of the wave vectors of 
emitted phonons must vanish (provided that the impu- 
rity scattering is neglected). For the same reason Raman 
scattering on one intervalley phonon must be impurity- 
assisted [process (b) in Fig. 1, giving rise to the so-called 
D Raman peak]. D peak is absent in the experimental 
Raman spectrum of graphene,^ showing that impurity 
scattering is indeed negligible in these samples. 

Looking at the intermediate electronic states involved 
in the Raman scattering (Fig. 1), we notice that for one- 
phonon scattering [processes (a), (b)] at least one inter- 
mediate state must be virtual, since energy and momen- 
tum conservation cannot be satisfied simultaneously in 
all processes. Thus, at least one of the factors in the de- 
nominator of Eq. (1) must be of the order of the phonon 
frequency ujph [for the impurity assisted scattering one 
of the electron-phonon matrix elements in the numerator 
of Eq. (1) should be replaced by the electron- impurity 
matrix element]. For the two-phonon scattering [pro- 
cess (c)] all intermediate states can be real, so that all 
energy mismatches in the denominator of Eq. (1) can be 
nullified simultaneously and the result is determined by 
the electron scattering rate 27. We emphasize the qual- 
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The full resonance picture presented above assumes 
the mirror symmetry between the electron and the hole 
spectra. The electron-hole asymmetry can be included 
as a correction to the Dirac spectrum: the electron and 
the hole energies can be written as vp ± aop^, where 
p is the momentum counted from the Dirac point, v is 
the Dirac velocity, and ao is the asymmetry param- 
eter. The energy scale Ae/t, quantifying the role of 
the asymmetry in the Raman scattering is defined as 



^eh = aoitof 



j)/(2?;)^, where Wi„ and uJout are the 



frequencies of the incident and the scattered photon (the 
details are given in Sees. VI B, VII C). Namely, the ar- 
guments of the previous paragraph hold if A^h <C 7. In 
the opposite case, it is A^h that determines the small- 
est value of the denominators in Eq. (1). We will always 
assume that both 7, A^h ^ (^ph- 

In the real space, the typical size of the region of space 
probed by the electron-hole pair in the fully resonant two- 
phonon Raman scattering, is ~ v/ max{7, Aeh}- For the 
doubly resonant defect-induced one-phonon scattering, 
the inverse energy mismatch l/tOph determines the time 
duration of the process by virtue of the uncertainty prin- 
ciple, so the length scale of the process in the real space 
is v/uiph- Although this length scale is much shorter than 
v/j, it is still much greater than the lattice constant or 
the electron wavelength v/e. 



B. Quasiclassical real-space picture 



FIG. 1: (color online). Schematic representation of the role of 
electron dispersion (Dirac cones, shown by solid lines) in the 
one-phonon (a,b) and two-phonon Raman scattering (c,d). 
Vertical solid arrows represent interband electronic transi- 
tions accompanied by photon absorption or emission (photon 
wave vector is neglected), dashed arrows represent phonon 
emission, the horisontal dotted arrow represents the impurity 
scattering. 



itative difference between the fully resonant process (c) 
and the double-resonant^^ process (b), where one inter- 
mediate state is still virtual. We also note the analogous 
difference between the two-phonon processes (c) and (d) 
in Fig. 1: only the process (c) is fully resonant, the other 
one involves an energy mismatch of 2LLjph- As a result, its 
amplitude will be smaller by a factor ^ j/LUph- 

Obviously, these arguments can be extended to all 
multi-phonon processes with odd and even number of 
phonons involved: in order to annihilate radiatively, the 
electron and the hole must have opposite momenta; if the 
total number of emitted phonons is odd, the electron and 
the hole must emit a different number of phonons, which 
is incompatible with energy conservation in all processes. 
Our main focus will thus be on even-phonon processes, 
as their intensities are determined by the electronic scat- 
tering. 



The fully resonant Raman scattering, where the en- 
ergy is conserved in each of the elementary scattering 
processes, admits a simple quasiclassical description, de- 
scribed qualitatively in this subsection, and justified rig- 
orously in Sees. VI, VII. 

Let us denote by e the energy of the electron and the 
hole in the photoexcited pair. Initially, it is given by the 
half of the excitation frequency Uin, e — ujin/2 ~ 1 eV. 
After the emission of ji phonons it is decreased by nojph] 
assuming Uph < 0.2 eV ^ e, we neglect this decrease in 
the qualitative considerations. Thus, during all the time 
taken by the Raman scattering, electron and hole can be 
viewed as wave packets of the size ^ w/e, propagating 
across the crystal along classical trajectories. 

The electron and the hole are created in the same re- 
gion of space of the size ~ v/e around some point ro at the 
moment of the arrival of the excitation photon. At this 
initial moment they have opposite momenta p, — p, and 
opposite velocities v, — v (if the electron-hole asymmetry 
is taken into account, the two velocities will have slightly 
different magnitude), so they move along the straight 
lines, their positions being Te{t) = ro+vt, r/i(i) — ro— vi. 
After a typical time t '^ I/7 they undergo some scattering 
processes (e. g., phonon emission), where their momenta 
and (generally speaking) energies are changed. Each such 
elementary scattering process occurs during a short time 
^ 1/e <C 1/7. Thus, the trajectories of the electron and 
the hole after their creation are represented by broken 



(a) 




(b) 
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FIG. 2: (color online), (a) An example of a quasiclassi- 
cal electron-hole trajectory contributing to the four-phonon 
Raman scattering. (b,c) Trajectories with emission of two 
phonons, not contributing (b) and contributing (c) to the 
two-phonon Raman scattering. In all pictures the lightning 
represents the incident photon which creates the pair. The 
solid lines denote the free propagation of the electron and 
the hole. The flash represents the radiative recombination of 
the electron-hole pair. The dashed lines denote the emitted 
phonons. 



II. SUMMARY OF THE MAIN RESULTS 
A. On the labelling of Raman peaks in graphene 

For the single-phonon Raman peaks the commonly ac- 
cepted notations are "G" for the peak at 1580cm~^ corre- 
sponding to emission of an optical phonon with zero wave 
vector, and "Z?" for the defect-induced peak at 1350cm~^ 
corresponding to emission of an optical phonon with the 
wave vector near K or K' points of the Brillouin zone.^ 
Sometimes one also distinguishes the so-called Z3'-peak 
at 1620 cm~^. This peak is also defect-induced, and cor- 
responds to emission of an optical phonon with a small 
wave vector q ^ e/v. As mentioned above, in the present 
work we study only the clean graphene, hence D and D' 
peaks are of no interest to us. 

Unfortunately, there is no single commonly accepted 
system for labelling of the multiphonon Raman peaks. 
The strong peak at 2700 cm~^, corresponding to emission 
of two phonons with the opposite wave vectors near the 
K and K' points, was historically called G" (as it is not 
defect- induced); sometimes it is denoted by 2D or by D* 
to stress that it is the second overtone of the D peak. 
The peak at 3250 cm~^ corresponding to emission of two 
phonons with two opposite wave vectors near the T point 
is sometimes called G*, 2G or 2D'. The latter notation 
reflects the fact that the frequency of this peak is not 
exactly the double of that of the G peak, but rather the 
double of the defect-induced D' peak. 

In the following we use the notation nT + mK to de- 
note the peak corresponding to emission of n phonons 
with wave vectors within ~ e/v from the F point and 
of m phonons with wave vectors within ~ e/v from the 
K or K' points. For multiphonon peaks this nomencla- 
ture is unambiguous. Thus, the peaks at 2700 cm"-'^ and 
at 3250 cm^^ will be called 2K and 2F, respectively. 



B. One-phonon Raman processes 



lines, with the typical segment length ~ v/^ (the elec- 
tron mean free path). The crucial point is that in order 
to recombine radiatively and contribute to Raman sig- 
nal, the electron and the hole should meet again within 
a spatial region of the size ^ t;/e, and have opposite mo- 
menta. The latter condition automatically implies that 
the number of the phonons emitted by the electron and 
the hole is the same. These considerations are illustrated 
by Fig. 2. 

In the presence of a significant electron-hole asymme- 
try, Ae/i :» 7, the described picture is modified. Namely, 
one of the segments of either the electron or the hole 
trajectory has the length w/Ae^ instead of w/7, the cor- 
responding time travelling being restricted by the phase 
mismatch rather than by collisions. 



In the clean graphene the only one-phonon Raman pro- 
cess, allowed by the momentum conservation corresponds 
to the emission of the E2 optical phonon with zero wave 
vector and frequency 1580 cm~^. For this process the 
situation turns out to be drastically different from that 
described by the qualitative considerations of Sec. I A. 
As shown in Sec. V, if one approximates the electron 
spectrum by the Dirac cones, the numerator of Eq. (1) 
vanishes due to high symmetry of the low-energy elec- 
tronic Dirac hamiltonian (as compared to the microscopic 
symmetry of the crystal). Thus, the main contribution 
to the Raman amplitude comes from the regions of the 
electronic Brillouin zone far from the Dirac points. As a 
consequence, the typical energy mismatch in the denom- 
inator of Eq. (1) is of the order of the whole electronic 
bandwidth. Thus, the Raman process, responsible for 
the 1580 cm~^ peak, is completely off-resonant and the 
picture shown in Fig. 1 for the process (a) is wrong. 



As a result, the intensity of the peak is expected to 
be insensitive to most external parameters: polarization, 
electron concentration, degree of disorder, etc. To char- 
acterize this intensity, one has to introduce an additional 
parameter into the theory which has no simple relation 
to the parameters of the low-energy effective hailtonian. 
The resulting intensity of the peak is given by Eq. (55); 
it is proportional to the fourth power of the excitation 
frequency which is the standard result for Raman scat- 
tering when the difference between the frequencies of the 
incident and scattered photons is small. This dependence 
also agrees with the experimental results of Ref. 14. 

Note that the results described above do not hold for 
the defect-induced peak at 1350 cm^^. For this peak the 
double resonance picture, ^^ shown in Fig. 1, process (b), 
is fully adequate. 



C. Two-phonon Raman processes 

As the phonons are emitted by electrons with momen- 
tum ujin/{2v), the largest possible phonon momentum is 
Qmax — i^in + ^out)/'2-v^ Corresponding to the electron 
and hole backscattering (for the 2K peak at 2700 cm~^ 
we count the phonon momenta from the K and K' 
points). It would be natural to expect that any pair of 
phonons with opposite momenta q, -— q and |q| < qmax 
can be emitted, the only exception being q = which 
is prohibited by symmetry^ and the nearby ones which 
are suppressed due to the smallness of the matrix el- 
ements. These arguments would predict the width of 
the peak to be of the order of {vph/v)u>in, where Vp^ is 
the K phonon group velocity; besides, the shape of the 
peak would be strongly asymmetric: a sharp cutoff on 
the high-energy side at the frequency 2ujpfi{qmax) due 
to the resonance restriction, and a smooth drop-off to- 
wards zero at 2ujph{q = 0) due to the matrix clement 
suppression. The phonon dispersion can be deduced 
from the dependence of the frequency of the impurity- 
assisted one-phonon D peak in graphite on the excita- 
tion energy ojin- dujph/duji„ = Vph/v ~ 50cm~^/eV.^^""'^® 
Thus, for LUin — 2 eV these arguments give the width of 
the 2K peak to be about 200 cm^^. However, the ex- 
perimentally observed width is only about 30 cm~^ at 
ujin = 2.2 eV, and its shape is quite symmetric. ^'^ 

The observed small width of the peak is explained by 
the quasiclassical picture, presented in Sec. IB. If upon 
the emission of phonons the electron and the hole are 
scattered by an arbitrary angle, as shown in Fig. 2(b), 
they will not be able to meet at the same spatial point in 
order to recombine radiatively and contribute to the two- 
phonon Raman peak. Only if the scattering is backwards, 
this event is possible, as illustrated by Fig. 2(c). This 
condition fixes the wave vectors of the emitted phonons 
to be g = qmax — 2e/v. The small deviations of the scat- 
tering angle from tt are restricted by the quantum diffrac- 
tion, and the width of the two-phonon Raman peaks, in- 
stead of being ^ {vph/v)ujin, is determined by a much 



smaller energy scale (see the discussion below). 

The dominance of the electron and hole backscattering 
manifests itself in the polariazion memory of the Raman 
signal. If the incident light is linearly polarized, the prob- 
ability of excitation of the electron-hole pair with a given 
direction of momenta is proportional to sin (p, where (p is 
the angle between the electric field vector of the light and 
the momenta. Thus, upon backscattering and radiative 
recombination, the probability to detect a photon of the 
same polarization as the original one is oc sin (p, and that 
of the orthogonal polarization is ex sin"^ (pcos^ (p. Aver- 
aging over (y9, we obtain the ratio of intensities for the 
detection of polarization parallel and perpendicular to 
that of the incident light to be I\\/I± = 3. This ratio 
may be slightly decreased due to a finite aperture (see 
the discussion in Sec. IV B). 

The calculation of the intensities of the two-phonon 
Raman peaks is performed in Sec. VI. The explicit ex- 
pressions for the intensities of the 2K and 2r peaks, ob- 
tained under the assumption of Dirac spectrum for the 
electrons, are represented by Eqs. (66) and (69). Both 
are proportional to 1/7^, where 27 is the electron (hole) 
inelastic scattering rate. If the latter is smaller that the 
electron-hole asymmetry A^h, then, according to the ar- 
guments of Sec. I A, it is Aeh that restricts the energy 
denominators from below. Formally, this results in the 
replacement (75) in both Eqs. (66), (69). 

Numerically, ao{l cVY/v^ ^ 0.1 eV (see, e. g., 
Ref. 19), so the relative correction to Eq. (66) for 
small Ae/i can be estimated as — (l/2)Ag^/(27)]^/2 ^ 
— 10~'*(a;j„/27)^. The total electronic broadening 27 was 
measured by time-resolved photoemission spectroscopy 
to be 20 meV in Ref. 20 and 25 meV in Ref. 21 (all values 
taken for e = oJin/2 — 1 eV). A recent angle-resolved pho- 
toemission spectroscopy (ARPES) measurement gives a 
significantly larger value for 27 ^ 100 meV (Ref. 22). 
Thus, the case 7 3> Aeh seems to be more relevant for 
the description of experiments, than the opposite one. 

The Raman matrix element corresponding to emission 
of two phonons with given wave vectors q and — q is given 
by Eq. (64) for 7 ^ Aeh- From this dependence one can 
deduce the lineshape of the two-phonon peaks 2K^ 2T: 
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(2) 



where /i = K,T ^ 2a;^ is the central frequency for each 
peak, and Vpu^^ is the group velocity of the corresponding 
phonon. Thus, the full width at half maximum (FWHM) 
of each peak is given by 



FWHlVhn = V2273 



1^^^^87; 



0.77^^2^87. (3) 



For 7 ^ Aeh the dependence of the Raman matrix ele- 
ment on q is described by Eq. (73). The lineshape cor- 
responds to two peaks separated by {vph,fj,/v)4:Aeh- Ex- 
perimentally, one sees just one 2K peak with the FWHM 
about 30 cm~^ at the excitation frequency ojin « 2eV.^'^ 
This corresponds to an unrealistically large value of 



27 ~ 0.2 eV. Most likely, this indicates that two-phonon 
peaks are broadened by other mechanisms, not taken into 
account in the present work. In particular, Eq. (2) ne- 
glects (i) the broadening of the phonon states, and (ii) the 
anisotropy of the phonon dispersion (trigonal warping of 
the phonon spectrum) . A detailed study of these effects 
would require introduction of additional parameters into 
the theory, so we prefer to postpone such study for the 
future work. It is worth emphasizing again that the in- 
tegrated intensity of the peaks, which is the main focus 
of the present study, does not depend on these details. 

In view of the results of the present paper it is worth 
mentioning the experimental measurements of the in- 
tensity I2K as a function of doping. While in Rcf. 9 
no significant dependence was observed, Ref. 12, where 
higher doping levels were reached, shows quite a strong 
dependence of I2K/ It on doping. The intensity /r of 
the off-resonant single-phonon 1580 cm^^ peak should 
not depend on doping (although the phonon width does 
exhibit such a dependence, the total spectral weight of 
the phonon state, determining the integrated intensity 
of the peak, must be preserved). At the same time, 
the intensity l2Ki if determined by the electron inelas- 
tic lifetime, should be sensitive to the concentration of 
carriers. Indeed, in the intrinsic graphene at low tem- 
peratures the photoexcited carriers do not participate in 
electron-electron collisions, as the phase space volume is 
restricted. ^^ As the carriers are added to the system, the 
electron-electron collisions become possible, thus the to- 
tal 7 increases, and the intensity I2K is decreased, in 
qualitative agreement with the observation of Ref. 12. 



D. Four-phonon Raman processes 

The motivation to study the four-phonon Raman pro- 
cess comes from the following picture for the fully reso- 
nant processes. The incident photon creates an electron 
and a hole - real quasiparticles which can participate 
in various scattering processes. If the electron emits a 
phonon with a momentum q, the hole emits a phonon 
with the momentum — q, and after that the electron and 
the hole recombine radiatively, the resulting photon will 
contribute to the two-phonon Raman peak. If they do 
not recombine at this stage, but each of them emits one 
more phonon, and they recombine afterwards, the result- 
ing photon will contribute to the four-phonon peak, etc. 
Three-phonon processes, not being fully resonant, are not 
interesting in this context. 

Besides phonon emission and radiative recombination, 
electron and hole are subject to other inelastic scattering 
processes, which can also be viewed as emission of some 
excitations of the system. In principle, Raman spec- 
trum should also contain the contrubution from these 
excitations, which are left in the system after the radia- 
tive recombination of the electron and the hole. The 
key point is that for real quasiparticles, the probability 
to undergo a scattering process a is determined by the 



ratio of corresponding scattering rate 27^ to the total 
scattering rate 27 = 'Ylia'^lai not by the history. This 
probability determines the relative frequency-integrated 
intensity of the corresponding feature in the Raman spec- 
trum. Thus, the ratio of integrated intensity I{2n+2)K of 
the Raman peak corresponding to 2n -\-2 K phonons to 
that for 2n K phonons {l2nK) niust be proportional to 
(7/^/7)^, where 2^k is the rate of emission of each of the 
two K phonons, and the square comes from the phonon 
emission by the electron and the hole. This conclusion 
depends weakly on the relation between 7 and Ae^ , only 
through a logarithmic factor. 

In the doped graphene, the most obvious competitor of 
the phonon emission is the electron-electron scattering: 
the optically excited electron can kick out another one 
from the Fermi sea, i. e., to emit another electron-hole 
pair. Thus, Raman spectrum should contain contribu- 
tion from electron-hole pairs; however, their spectrum 
extends all the way to the energy of the photo-excited 
electron (optical energy) in a completely featureless way. 
Thus, it cannot be distinguished from the parasitic back- 
ground which is always subtracted in the analysis of Ra- 
man spectra, and cannot be seen in the Raman spectrum 
directly. However, assuming 7 = ^k + 7r + 7ee , where 
27r is the rate of emission of phonons from the vicinity 
of the r point of the first Brillouin zone, and 27ee is the 
electron-electron collision rate, one can extract the value 
of 7ee, relative to phonon emission rates from the experi- 
mental data. More precisely, in this way one obtains the 
rate of all inelastic scattering processes where the elec- 
tron loses energy far exceeding the phonon energy. 

Note that arguments leading to I(2n+2)K/l2nK « 
{iphllY are not specific for graphene; in fact, this is 
nothing but Breit-Wigner formula, applied once for the 
electron and once for the hole. Multi-phonon Raman 
scattering has been studied in wide-gap semiconductors 
both experimentally^^'^^ (up to ten phonons were seen 
in the Raman spectra of CdS), and theoretically^^'^^. In 
a wide-gap semiconductor an optically excited electron 
does not have a sufficient energy to excite another elec- 
tron across the gap, so the electron-electron channel is 
absent. In addition, interaction with only one phonon 
mode is dominant, so the ratios of subsequent peaks are 
represented by a sequence of fixed numbers. The sim- 
ple band structure (one valley for CdS in contrast to two 
valleys for graphene) allowed a calculation of the whole 
sequence. A more complicated electronic band structure 
in graphene makes it problematic to calculate the whole 
sequence, so wc restrict ourselves to the calculation of I^^k 
for the most intense four-phonon peak. 

This calculation is performed in Sec. VII. Its result 
depends, besides the relation between 7 and Ag/i, also 
on their relation to the energy scale u>in{vph/v), char- 
acterizing the phonon dispersion. In Sec. II C we have 
already discussed this energy scale; for Wm = 2 eV we 
have uJin{vph/v) « 100 cm~^ « 12 meV. The meaning of 
this energy scale is the difference between the energies of 
the electron and the hole after each of them has emitted 



two phonons with almost arbitrary momenta (the only 
restriction is that the sum of all four phonon momenta 
must vanish). If ujin{vph/v) <C 7, A^h, which seems to be 
the case (see the discussion in the previous subsection) 
then this difference can be neglected, and the intensity 
of the 4K peak is given by Eq. (101) for 7 3> A^h [which 
is likely to be the case relevant for most experiments, 
and which was reported in the short paper by the author 
(Ref. 28)], and by Eq. (104) for 7 < A^h; the polariza- 
tion memory is lost in both these cases, /|| « J^. In 
the case tUinivph/v) 3> j,Aeh the intensity I4K is given 
by Eq. (118), and a significant polarization memory is 
expected, up to /|| /I± « 3. 

A thorough experimental study of the intensity I/^k 
(in particular, its dependence on doping) is still lacking. 
Our prediction for the case 7 ^ A(,}i,ujin{vph/v), which 
we believe to be the experimentally relevant one, is^*'^' 



I2K 



0.11 



iK 



"fK +"fr + "fee 



(4) 



E. Renormalization of the coupling constants 

In the calculations of the Raman intensities, described 
above, electron-phonon coupling constants entered as pa- 
rameters of the theory, without any assumptions about 
their values, except for the relations fixed by the sym- 
metry of the crystal. In particular, the two-phonon peak 
intensities I2K and /2r are determined by two indepen- 
dent dimcnsionless coupling constants which we denote 
Xk and Ar [see Eq. (24) for the definition]. A simple 
estimate of the coupling constants can be obtained from 
the tight-binding nearest-neighbor model of the graphenc 
crystal. In this model the only parameter characterizing 
the electron spectrum is the nearest-neighbor electronic 
matrix element to, and the electron-phonon interaction is 
characterized by its change with the bond length, dto/da. 
In this model we obtain Xk /^r to be given by the inverse 
ratio of the corresponding phonon frequencies, about 1.2; 
the same result up to a few percent is obtained from the 
density-functional theory (DFT) calculations of Ref. 29. 
At the same time, by comparing the experimentally mea- 
sured intensities of the different two-phonon peaks, and 
using the result of our calculation performed in Sec. VI, 
we can independently extract the ratio of the coupling 
constants. According to the data of Ref. 6, I2K /I2V ~ 20, 
which gives Xk /Xy ~ 3. 

To explain this discrepancy we first analyzed the ef- 
fect of the electronic trigonal band warping, which affects 
I2K and /2r differently. The corresponding calculation is 
done in Sec. VI B. For ujin = 2 eV, we estimate the rela- 
tive contribution of the warping term as 5 • 10~ for I2K 
and 5 • 10^^ for l2r, which is far too little to account 
for the observed ratio hK/hr- We are thus led to the 
conclusion that the observed ratio I2K/I2T must be due 
to the difference of the coupling constants, not accounted 
for by the DFT calculation. A similar conclusion about 



the insufficiency of the DFT calculation of the electron- 
phonon coupling constants has been drawn in Ref. 31, 
where an attempt was made to explain the experimental 
data obtained by ARPES^^ using the results of the DFT 
calculation. 

At the same time, we should note that the dimcn- 
sionless coupling constant Ar for the phonons near the 
F point, as calculated by DFT^^ (Ar ~ 0.028), agrees rea- 
sonably well with the measured one: the measurements 
of the linear in the wave vector q term in the phonon 
dispersion (Kohn anomaly due to electron-phonon inter- 
action), ujph(q) - ujph{q == 0) w {Xr/8)vq, give Ar w 0.024 
(see Ref. 30) ; the measurements of the dependence of the 
phonon frequency LUph on the electron Fermi energy ep' 
Aujph ~ (Ar/27r)|eF| give Ar w 0.034 (Ref. 9) Ar w 0.027 
(Ref. 10). 

We show that the difference between the ratio 
^K / ^r ~ 3 extracted from the Raman peak intensities, 
and Xk/Xy w 1.2 obtained by the DFT calculation,^^ is 
due to the part of Coulomb interaction between electrons, 
not picked up by the DFT when local approximations 
are used for the exchange-correlation functional, such as 
the local density approximation (LDA) or the general- 
ized gradient approximation (GGA), namely, logarith- 
mic renormalizations.'^^ Coulomb interaction has been 
known to be a source of logarithmic renormalizations 
for Dirac fermions.'^^""^^ Coulomb renormalizations in 
graphene subject to a magnetic field have been consid- 
ered in Ref. 35, Coulomb effect on static disorder has 
been studied in Refs. 36-38. Essentially, the idea of the 
renormalization of the coupling constants is that the ma- 
trix element of the electron-phonon interaction should be 
taken not between the non-interacting electronic states, 
but between the states dressed by the Coulomb interac- 
tion. If the typical electronic energy in the problem is e 
(~ 1 eV in the case of Raman scattering) , the renormal- 
ization is determined by the Coulomb interaction at all 
length scales from the shortest ones (lattice constant) to 
the electron wavelength v /e. It is this long-range part 
of the exchange and correlation that is missed by the lo- 
cal approximations in the DFT calculation, which take 
into account correctly only the short-range correlations 
(at the distances of the order of the lattice constant). 

In Sec. VIII A we calculate the renormalization of the 
dimcnsionless electron-phonon coupling constants (a pre- 
liminary account of this work was given in the short 
paper'^^), and show that the coupling constant Ar for the 
phonons near the F point is not renormalized (hence the 
agreement between the value of Ar calculated by the DFT 
and measured in the experiments, as mentioned above), 
while the coupling constant Xk for the phonons near the 
K point, which is responsible for the 2K Raman peak, is 
enhanced by the Coulomb interaction. This enhancement 
depends on the electronic energy, as shown in Fig. 22. For 
the electronic energy 1 eV this enhancement is in quanti- 
tative agreement with the measured ratio I2k/I2Ti pro- 
vided that the screening of the Coulomb interaction by 
the substrate is weak. The dependence of the enhance- 



ment on the electronic energy translates into the depen- 
dence of I2K I I2T on the excitation frequency, which can 
be checked experimentally. Similarly, as the Coulomb in- 
teraction is screened by the substrate with a dielectric 
constant £00 (its high-frequency value), the dependence 
of I2K / I2T on £00 can also serve as an experimental check 
of the theory. 

We also show in Sec. VIII B that the electron-phonon 
coupling itself is a source of logarithmic renormalizations. 
However, due to the smallness of the coupling constants 
this effect is much weaker than the effect of the Coulomb 
interaction. 



F. Structure of the paper 



In Sec. Ill the low-energy hamiltonian of the interac- 
tion of electrons with the crystal vibrations and the elec- 
tromagnetic field is written from pure symmetry consid- 
erations. In Sec. Ill A the symmetry of the graphenc 
crystal is reviewed. In Sec. IIIB the symmetry con- 
siderations are used to write the electronic part of the 
hamiltonian. Sec. IIIC is dedicated to the symmetry 
analysis of the Dirac part of the electron hamiltonian, 
whose symmetry is significantly higher than the symme- 
try of the crystal. Sec. Ill D is dedicated to the symmetry 
analysis of the in-plane crystal vibrations. In Sees. HIE 
and III F we write the hamiltonian of interaction of elec- 
trons with the optical and acoustical vibrations, respec- 
tively. Sec. IIIG is dedicated to the symmetry analysis 
of the out-of-plane vibrations of the graphene crystal. In 
Sec. Ill H the hamiltonian of the interaction of electrons 
with the electromagnetic field is written. 

Sec. IV describes the general scheme of the calculation 
of Raman peak intensities using the standard perturba- 
tion theory. In Sec. IV A the Green's functions are in- 
troduced, and in Sec. IV B the general expression for the 
Raman scattering probability is derived. In Sec. IV C 
we discuss the electron inelastic scattering, and calcu- 
late the electronic self-energy due to the electron-phonon 
coupling. 

In Sec. V one-phonon Raman scattering is discussed, 
and it is shown that the calculation of the one-phonon 
peak intensity cannot be performed within the low- 
energy theory. In Sec. VI the two-phonon Raman peak 
intensities are calculated, first, under the assumption of 
the Dirac electron spectrum (Sec. VIA), and then tak- 
ing into account the trigonal band warping and electron- 
hole asymmetry (Sec. VI B). In Sec. VII the intensity 
liK of the most intense four-phonon peak is calculated. 
Sees. VIII A, VIII B are dedicated to the renormaliza- 
tion of the electron-phonon coupling constants due to 
Coulomb interaction and due to the electron-phonon in- 
teraction, respectively. 




FIG. 3: The honeycomb lattice with two atoms [A and B) per 
unit cell. Tripled unit cells are shown by dashed hexagons. 
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TABLE I: Irreducible representations of the groups Ceu and 
Csv and their characters. 



III. SYMMETRIES AND HAMILTONIAN^" 
A. Symmetry of the crystal 

Since the typical energy of the incident photon (about 
2 eV) is much smaller than the 7r-electron bandwidth 
(~ 20 eV), one can expect the low-energy excitations 
to play the dominant role. In this section we employ 
standard symmetry analysis^^ to fix the form of the low- 
energy hamiltonian. We prefer not to choose any specific 
basis and use algebraic properties. 

The carbon atoms of graphenc form a honeycomb lat- 
tice with two atoms per unit cell, labeled A and B 
(Fig. 3), the distance between nearest neighbors being 
a = 1.42 A. Three out of four electrons of the outer shell 
of each carbon atom form strong a bonds with its three 
nearest neighbors, and represent no interest to us. The 
remaining tt orbitals (one per each carbon atom) give rise 
to the half-filled tt band. 

In this paper we will not consider the dimension, per- 
pendicular to the crystal plane, so the point symmetry 




FIG. 4: The first Brillouin zone corresponding to the honey- 
comb lattice, and its wrapping upon triphng of the unit ceU. 
Regions with the same shading should be translated so that 
the K, K' points move to the F point to form the first Bril- 
louin zone of the crystal with the tripled unit cell. Solid lines 
show the (Ta,ai,,(Jc refiection axes. The a„,CT^,a^ refiection 
axes are shown by dashed lines (not labeled). 



group of the graphcne crystal is C'sy. It contains 12 
elements: the identity, five rotations Cq, n = 1,...,5 
{Cm denoting the rotation by I'k jm) and six reflections 
in planes perpendicular to the crystal plane. The three 
reflections leaving the A and B sublattices invariant are 
denoted by da, CTf,, CTc, while those svifapping the A and B 
sublattices points will be denoted by u'^^^a'^^^a'^. Table I 
lists the irreducible representations of the group Cgt; and 
their characters. 

The flrst Brillouin zone of the crystal is a hexagon 
(Fig. 4). Out of the six corners of the hexagon only two 
are inequivalent. They are called K and K' points. The 
group of the wave vector at these points is Czv- The 
states at these points are twice degenerate, transforming 
according to the two-dimensional irreducible representa- 
tion (i?) of C-iy. Transformations, swapping K and K' ^ 
thus belonging to Cgi, but not to Czv , (reflection a^ in the 
plane, perpendicular to that of cr(,, and rotations C2, Cg), 
flx the energies at K^ K' to be equal. One can form real 
linear combinations of wave functions from K and K' 
points which transform according to Ei and E2 repre- 
sentations of Cqv The degeneracy at K,K' points, in 
combination with the absence of any other states in the 
Brillouin zone with the same energy, flxes the Fermi level 
of a half-flUed band to be at this energy, which is thus 
natural to choose as e = 0. 

Instead of dealing with degenerate states at two differ- 
ent points of the Brillouin zone {K, K') , one can triple 
the unit cell of the crystal. The new unit cell contains 6 
atoms which form a hexagon (Fig. 3), while the new Bril- 
louin zone is only 1/3 of the original one (Fig. 4). The 
advantage of this approach is that now both K and K' 
are mapped onto the F point, so one does not have to 
consider the two of them separately. 

Tripling of the unit cell means that two translations 
tai and ia2 are factorized out from the translation group 



of the crystal, so they should be added to the point group, 
which becomes Cg'y = Cqv + {ta^CQ^) + {ts^^Cev). Irre- 
ducible representations of this group and their charac- 
ters are shown on Table II. The states of the 7r-electrons 
at the new F point form a 6-dimensional representation 
which is reduced as A1+B2 + G' , where Ai and B2 states 
are the non-degenerate ones corresponding to the old 
F point, while the 4-dimensional irreducible representa- 
tion G" contains the zero-energy states inherited from the 
old K, K' points. 

In order to write down the low-energy electronic hamil- 
tonian, we have to consider 4x4 hermitian matrices act- 
ing in the 4-dimensional space of the zero-energy elec- 
tronic states. The basis in the 1 6-dimensional space of 
such matrices is provided by the generators of the S'C/(4) 
group forming a 16-dimensional reducible representation 
of Cqv or Cg'y . This representation is reduced as 

{Ei+E2)x{Ei + E2) = 2iAi + A2 + Bi + B2 + Ei + E2), 

(5) 
within the group Cgt, (the two sectors corresponding to 
matrices either diagonal or off-diagonal in the KK' sub- 
space), or as 

G'xG" == Ai+Bi+B2+A2+Ei+E2+E[+E'^ + G', (6) 

within the group Cq^. The correspondence between 
Eqs. (5) and (6) is given in Table III. 

The most convenient way to identify the matrices is 
by specifying the irreducible representation according to 
which they transform, rather by specifying their explicit 
form in some particular basis. So, the matrix which 
transforms according to the A2 representation of Cg'y 
will be denoted by S^ and called the z-component of 
the isospin (the only arbitrariness in this deflnition is the 
overall sign). The two matrices which transform accord- 
ing to the vector Ei representation of Cg^ will be denoted 
by {Sj;, Ey} = S [defined up to an arbitrary rotation, see 
Eq. (15a) below], and so on. The full list of definitions 
and notations is given in Table III. 

Explicit expressions for the electronic matrices are not 
needed, as long as their algebraic rules are specified. 
The simplest way to specify these rules is to express 
all the 16 matrices in terms of two sets, {Sa;,Sj,,Sz} 
and {Ax, A.y,Az}, and their products, where the matri- 
ces from the same set satisfy the Pauli matrix algebra, 
while matrices from different sets just commute. In Ap- 
pendix A we show how these rules can be established, and 
give the explicit expressions for the matrices for some spe- 
ciflc choices of the basis. We also note that q(^-^^/^)^^ js 
the matrix of the C3 rotation, A^'S^ - of the C2 rotation, 
AzSa; - of the CTq reflection, AyY,y - of the CTq reflection, 
and e^i'^'^^/'i)^^ are the matrices of the two elementary 
translations. 



B. Electronic wave functions and hamiltonian 

Since there are two 7r-orbitals per unit cell, the eigen- 
functions for each wave vector k in the Brillouin zone are 
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TABLE II: Irreducible representations of the group Cg^ = Cav + (taiCeu) + (ia2C'6u) and their characters. 
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TABLE III: Classification of 4 x 4 hermitian matrices and their transformation properties under time reversal T. 



given by e Uk_±{r, z), where the two Bloch functions 
^k.±(r, z), periodic in r = (x, y), are even and odd with 
respect to the reflection a'^^, and are normalized as 



jKAr,z)UUr,^)d^rdz^L,.L,S.,, 



^,] 



(7) 

where LxLy is the crystal area. It is more convenient 
to choose their linear combinations, localized near each 
carbon atom: Z^k,A — (^k,+ + ^k,-)/'\/2 and Z^k.s = 
(^k.+ — ^k.-)/v2- In this basis an arbitrary wave func- 
tion ^(r, z) involving only low-energy states can be writ- 
ten in terms of a four-component smooth envelope func- 
tion V'(r), or its Fourier transform ^'(p); pa <C 1, as 



^(y,z) = 



rf2p 



(27r 



+ V'K',A(p)e^<^'+P)"-Z^K',A(r,z) 
+ V^K'.B(p)e*<^'+P^'"WK',ij(r,z) . (8) 



The low-energy effective electronic hamiltonian i7oi(p) is 
defined as a 4 x 4 matrix whose matrix element between 
any two smooth envelope functions ^{y) and V'(r) coin- 
cides with the matrix element of the microscopic hamil- 
tonian 7iei(— *V, —idz', r, z) including the periodic crystal 
potential, between the corresponding full wave functions 



^(r, z) and ^(r, z) (see Appendix B for details): 

V'f(r)iJei(-iV)^(r)d2r = 



= / **(r,z)7i:ci(-iV,-i9^;r,z)*(r,z)d2rdz.(9) 



Strictly speaking, the spin index should also be at- 
tached to the enevlope function V'(r). However, it would 
make the formulas more cumbersome, and we prefer to 
omit it, as none of the calculations of the present paper 
will concern a non-trivial spin structure. The only role 
of the spin will be to provide an additional degeneracy, 
which will be accounted for, and mentioned separately 
every time it will enter the calculations. 

We expand the effective hamiltonian in powers of p: 
i7ci(p) - Hi{p) + H2ip) + . . ., where il„(p) = 0(p"). 
One can write down different terms from symmetry con- 
siderations, taking into account that momentum compo- 
nents PxiPy transform according to Ei (vector) repre- 
sentation of Cgi,. The leading term in the hamiltonian, 
-ffi(p), must have the Dirac form:"'^ 



i7i(p) =: wpj;Ej; +wpj,Ej, = wpE. (10) 

The coefficient v turns out to be equal to u w lO^cm/s w 
7 eV • A; it can be related to the nearest-neighbor cou- 
pHng matrix element t of the tight-binding model as 
V — 3ia/2. The four eigenstates of the hamiltonian (10) 
for each p can be classified by the value (±1) of projec- 
tion of the isospin S on p, corresponding to the energies 
±w|p|. In other words, the hamiltonian (10) is diagonal- 
ized by a unitary transformation: 
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(Pp — arctan(pj^/pa;) being the polar angle of the vector p. 

Various perturbations of the Dirac hamiltonian (10) 
should also be classified according to Table III. Only 
those containing the matrix S^ will open a gap in the 
electron spectrum. Perturbations, not containing Si, cor- 
respond to an energy shift of the whole spectrum (which 
may be accompanied by valley mixing, if A matrices are 
involved). Perturbations, proportional to S, correspond 
to a momentum shift of the Dirac points, which can be 
viewed as a gauge vector potential. 

The next term in the hamiltonian can be written by 
taking into account that the symmetric tensor pipj can 
be decomposed into p^, transforming according to Ai, 
and 2pxPy,pl — Py, transforming according to E2: 

H2{p) = ao/l + "3 [-2p^pyT.y + {pI-pI)T,^] a,. 

The first term describes electron-hole asymmetry; it van- 
ishes in the nearest-neighbor tight-binding model, and 
appears only if coupling to the second nearest neighbors 
is included. The second term, corresponding to E2 rep- 
resentation, is responsible for the so-called trigonal band 
warping. In the nearest-neighbor tight-binding model its 
value is given by — 3ta^/8. 

At this point it is convenient to introduce the time- 
reversal operation whose action on the microscopic spin- 
less wave function is defined by ^(r) 1-^ ^*(r). For 
the four-component envelope function ip(r) this defini- 
tion translates into 



symmetries are (i) the full intravalley rotational symme- 
try Coot,: 



V'(r)h-.C/rV'*(r), 



r 



(13) 



where Ut is a unitary 4x4 matrix, with an additional 
requirement Uq-U^ ~ 1, whose explicit form depends on 
the choice of the basis. The action of time reversal on 
the effective electronic hamiltonian is defined by 



H,i{p,r)^UTH:,{-p,r)ul 



T- 



(14) 



Behavior of electronic matrices under the time reversal 
is listed in Table III. 

Time reversal symmetry of the hamiltonian does not 
add any new symmetries to the spectrum, as compared to 
those imposed by the spatial symmetry Cqv (namely, C^v 
symmetry of the spectrum around each of K, K' points, 
and the mirror symmetry between the spectra at K and 
K' points), because the action of the time reversal on 
the wave vector is identical to that of the C2 rotation, 
K + p H^ K' — p. However, some perturbations may lift 
the C2 symmetry, while still preserving the time reversal 
one (see Sees. Ill E and III H) 



C. Additional symmetries of the Dirac hamiltonian 

The Dirac hamiltonian (10) has a higher symmetry 
than the microscopic symmetry Cq^. The additional 




-iE^cp/2-c-i iSjcp/2 

COS if — sin if 
sin if cos if 




(15a) 
, (15b) 



which leaves i?i(p) invariant; (ii) the chiral property: 

UcHiUc - -Hi, Uc = S„ (16) 

which ensures the symmetry of the spectrum with respect 
to e — > —e (i. e., particle-hole symmetry). 

The intravalley "time reversal" symmetry, mentioned 
in Ref. 43, can be represented as a combination of the 
time reversal (13), the intravalley rotation (15a) by tt, 
and the C2 rotation: 

^ ^ A,E,e-*^-"/2[/r^* = -iA^Ut^- (17) 

Applying this operation twice results in a minus sign, 
since the matrix A^; is odd under time reversal. 

If one wishes to include the second-order hamilto- 
nian (12), the Ai term preserves only the rotational sym- 
metry Cocv, while the E2 term preserves only the chiral 
property (16). 



D. In-plane phonon modes 

There is quite extensive literature dedicated to the 
phonon modes of graphene and graphite, and their sym- 
metry analysis (see, e. g., Ref. 2 and references therein). 
To make the presentation self-contained, we briefly re- 
peat the facts which are necessary for the subsequent 
considerations. 

The only phonons that can efficiently couple to the 
low-energy electronic states are those near the P point 
(coupling electronic states in the same valley), and near 
the K, K' points (coupling electronic states in different 
valleys; note that K — K' is equivalent to K'). As a re- 
sult, each unit cell has four in-plane degrees of freedom 
(two per each carbon atom). At P point one has two 
acoustic and two optical modes. Coupling of the acous- 
tical modes with a wave vector q to the electron motion 
must vanish as (7 ^ 0, and thus be small in the parame- 
ter qa, but we consider them for the sake of completeness. 
For eight modes at K, K' points it is more convenient to 
consider their real linear combinations which transform 
according to Ax + Bi + A2 + B2 + Ei + E2 representa- 
tions of the group Cgi,. All twelve modes are shown in 
Fig. 5. They are linearly polarized, in contrast to the 
basis with a definite wave vector, which would have a 
circular polarization. 

After having mixed the K and K' modes, it is nat- 
ural to switch to the tripled unit cell representation. 
Analogously to Eq. (8) for electrons, an arbitrary lat- 
tice displacement pattern, involving phonon states with 
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FIG. 5; In-plane phonon modes at F point and real linear combinations of the phonon modes at K, K' points corresponding to 
different irreducible representations of Csv The dashed lines show tripled unit cells. 
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wave vectors close either to F, or K, or K' points, can 
be expressed in terms of a smooth envelope function 
u{v) or its Fourier transform u(q). Since we have to 
specify two cartesian components of displacements for 
each of the six atoms in the tripled unit cell, 5x_ = 
((5xi, (^2/1, . . . , 5xq, Si/q)'^ , the envelope of normal coordi- 
nates u{q) is a twelve-component vector: 



fe(R)= /'^f;n^(q)e^<i^A:^ 



(18) 



where R spans a discrete set of points, labeling different 
unit cells, and X_ is the pattern of displacements for 
the /ith normal mode, normalized as 



X 



^q^prilq,^' = 6(5 



'^J.^J.'■ 



(19) 



For q = the twelve normal modes X n ^ are shown in 
Fig. 5. Introducing 11^ (r), the canonically conjugate mo- 
mentum to u^{v), we write the bare phonon hamiltonian 
as 



u— 1 ^ 



mv) , M 



2M 



^"a'(i")'^»(-*^)"mW 



(20) 
with M ~ 1.993 -10 ^'^ g being the carbon atom mass, 
N the total number of the carbon atoms in the crystal, 
and Wp(q) the frequency of the /ith normal mode. 

At q = the modes belonging to the same irreducible 
representation of Cg^ are degenerate. Expansion of the 

I 



phonon potential energy at small q, describing the split- 
ting, can be done analogously to that of the electronic 
hamiltonian. It is convenient to preserve the form (20) 
of the hamiltonian, but instead of summing over normal 
modes /i one should sum over irreducible representations 
of Cg'y. For each irreducible representation cj^(— «V) be- 
comes a matrix acting in the space of the degenerate 
modes belonging to this representation. 

For two-dimensional representations {Ei and E2 at the 
r point, and E'^^E'2 at K,K' points) we have to classify 
2x2 hermitian matrices according to irreducible repre- 
sentations in the corresponding decomposition: 



Eix Ei= A1+A2 

E2X E2= A1+A2 

E[xE[=Ai+A2 



E2, 


(21a) 


E2, 


(21b) 


E[, 


(21c) 


E'2- 


(21d) 



E'2xE'2 = Ai+ A2 



The basis in the space of such matrices is provided by 
the unit matrix ctq^ transforming according to Ai, the 
Pauli matrix a^ transforming according to A2, and the 
Pauli matrices Ox , (Jy transforming according to the third 
term in each decomposition. Since the components qx, Qy 
transform according to Ei , the linear in q term is absent 
for all two-dimensional representations. For the four- 
dimensional representation G" we have to deal with 4x4 
matrices, so everything is fully analogous to the electronic 
case, and these phonons also have Dirac spectrum. In the 
second order we have q'^ ^ A\ and q^ ~ Q^^ '^Qxqy ^ E2, 
so we can write 



wL(q) 



vl + vl 2 



[(9x - (il)'^x - 2qxqy(Ty 



= ^E2 + ^e2{Jl + JT)q^ + i^e2{Jl - Jt) [{ql - ql)<Jx - iq^qyUy] , 



w|;(q) 

w^,(q) = w^, +2wG'WG'qS + 0(g^ 



4; +2ujE'jE[q^, 
j%, +2ujE'JE!y, 



(22a) 

(22b) 
(22c) 
(22d) 
(22e) 



The 2x2 matrices appearing in Eqs. (22a), (22b), can 
be diagonalized to yield the dispersion of the longitu- 
dinal and transverse phonons: '-^Ei,L{T)i^) — vl{T)<1i 

WB2,i(T)(q) = ^E2 + JL(T)q^- 



E. Electronic coupling to in-plane optical phonons 

In the electron-phonon interaction hamiltonian the 
normal mode displacements should be paired with the 
electronic matrices, corresponding to the same irre- 
ducible representation. For optical phonons it is suffi- 
cient to take the leading q = term, so we have one 



independent coupling constant for each irreducible rep- 
resentation of the group C'^^ : 



iJo-opt = FT[nE2 X A^S]^ + 

-I- FK{uA^^-x^z+UB^l^y'i^z) + 

+ F'k {[vle, X A^S]^ + U£,Aj,S) 



(23) 



Note that E'2 phonons cannot couple to the electron mo- 
tion in this approximation since the corresponding matri- 
ces Aa;, Ay change sign under the time reversal, while the 
electronic hamiltonian must preserve time reversal sym- 
metry even in the crystal with displaced atoms. Electrons 
will couple to (i) momentum II^;' , whose effect is small 
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in the ratio of the phonon frequency u>e' to the electron 
bandwidth, (ii) gradient Vue' , whose effect is small in 
the parameter qa, (iii) the squares of displacements u'^^, , 

whose effect is small in the parameter l/{Ma'^ujE' )• 

The three constants in Eq. (23) can be evaluated in the 
tight-binding approximation, where they are expressed 
in terms of F = dt/da ss 6 eV/A - the change of the 
nearest-neighbor coupling matrix element with the dis- 
tance between the atoms. We obtain Fy =- Fk = 3F, 
while F'j^ vanishes. This vanishing is an artifact of 
the nearest-neighbor bond-stretching approximation. If 
one takes into account the dependence of the electronic 
hamiltonian on the angles between the bonds, F^ be- 
comes different from zero. Nevertheless, it is an order 
of magnitude smaller than Fk (e. g., density- functional 
theory calculations^^ give F^/Fk ~ 0.13), and we ne- 
glect the G" phonons for the rest of the paper. Thus, 
our attention will be focused on the two degenerate 
modes E2 and E[, which will be referred to simply as 
r and K phonons, and their frequencies will be de- 
noted by uJE2,L{q) == wr,L(q), w^^^tIq) = tJr,T(q), 
ujE[{q) = WK(q). 

The strength of elcctron-phonon interaction can be 
conveniently characterized by the dimcnsionlcss coupling 
constants 



Ar = 



F? 



27a' 



Mwr(0)u2 



Xk = 



P2 



27a' 



Mwk(0)w2 



(24) 

where \/27a^/4 is the area per carbon atom. For Fp — 
Fk = 3F, F = 6eV/A, M = 2.00-10"23g ^ 2.88-103(cV- 

A^)-i, V = 10*3 m/s = 6.58eV- A, wr = 1580 cm^^ = 
0.196 eV, UJK = 1370 cm-i = 0.170 eV, a = 1.42 A we 
have 



Fl V27a2 ^ 


« 0.035, 


Fl V27a2 


Mloyv' 4 


Mlokv' 4 



0.040. 

(25) 

The easiest way to match the notations for the electron- 
phonon coupling constants, used in the present paper, to 
those used in other works, is to compare observable quan- 
tities. For example, adding electrons to the system leads 
to a shift of the phonon frequencies. For the F phonons 
at q = in the present notations this shift is expressed 
in terms of the Fermi energy tp as 






eF 



cjr , loy-2\(.f\ 

—r In 77x — ; 

4 wr - 2 ei? 



(26) 



Alternatively, one can look at the correction to the 
phonon dispersion as a function of the wave vector q 
(see Appendix E). 



only through spatial and time derivatives of the corre- 
sponding displacements: 



n, 



M 



- (p + TTlef S) 



-f 'E.Q{dxUEix + dyUEiy)t + 

+ 'E.i{dxUE^y ~ dyUE^x)'^z + 

+ S2 [-{dxUE^y + dyUEix)^y + 

+ {dxUExx - dyUEj_y)^x] A^. (27) 

The general form of the coupling to TIei in first line fol- 
lows from the symmetry considerations (the coefficient at 
n^;^ should transform according to Ei and change sign 
under time reversal). However, its exact form follows 
from the k • p perturbation theory (which identifies Wg 
with the free electron mass) in combination with the re- 
quirement of Galilean invariance: if ^(r,f) satisfied the 
Schrodinger equation for the stationary lattice, then for 
the lattice moving at a constant velocity u^^ = Uei/M, 
the solution should be ^(r — li^j^i, t) (we neglect the con- 
tribution of electrons to the total kinetic energy of the 
crystal). Coupling to the strain is written using the de- 
composition El X El — Ai -\- A2 + E2- 



G. Out-of-plane phonons 

The basis vectors for out-of plane displacements coin- 
cide with those for the electronic wave function ampli- 
tudes in the tight-binding picture (since in both cases a 
single number is associated with each lattice site). At 
the F point the basis vectors represent A and B atoms 
shifting in the same direction or opposite directions, cor- 
responding to Ai and B2 representations. For the K, K' 
phonons we have Ei + E2 representations of Cq^, , or the 
four-dimensional G" representation of Cg^ in the tripled 
unit cell picture. All six modes are shown in Fig. 6. 

For a suspended graphene sheet the Ai mode has a fre- 
quency lliai{c\) oc q^ (see, e. g., Ref. 44). The frequencies 
of other modes are finite. All frequencies can be strongly 
affected by the interaction with the substrate (in par- 
ticular, there will be no reason for the frequency of the 
Ai mode to vanish at q = 0). 

Out-of-plane displacements can couple to the electron 
motion only quadratically, if the crystal is symmetric 
with respect to reflection in the crystal plane. While 
it is the case for a suspended graphene sheet, presence of 
a substrate may break this symmetry.^^ 



H. Coupling to electromagnetic field 



F. Electronic coupling to in-plane acoustical 
phonons 

Since a uniform translation of the crystal cannot affect 
the electron motion, it can couple to acoustic phonons 



The hamiltonian of interaction of electrons with the 
electromagnetic field, described by the long-wavelength 
scalar and vector potentials (^{y) and A(r), can be 
obtained from the requirement of gauge invariance: 
-^oi(p) -^ Hc\(p — (e/c)A) + e(p. To the linear order 
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FIG. 6: Out-of-plane phonon modes at F point and real lin- 
ear combinations of the phonon modes at K, K' points corre- 
sponding to different irreducible representations of Cev- 



in A we have 



Hr 



ev . .„ 

~ eif AS. 

c 



(28) 



Besides coupling via gauge potentials, electrons can 
couple directly to electric and magnetic field. Such terms 
are gauge invariant and cannot be deduced from the bare 
electronic hamiltonian in the envelope function approx- 
imation, as they correspond to the effect of electromag- 
netic field on the microscopic Bloch functions. They 
should be introduced into the effective low-energy the- 
ory directly, and can be either calculated microscopically 
(see Appendix B for the calculation by k • p perturbation 
theory), or written from symmetry considerations. 

The electric field vector is invariant under time re- 
versal, and its cartesian components £x,£y,£z transform 
according to £z ~ Ai, {£x,£y) ~ Ei under C^v The 
magnetic field vector changes sign under time reversal, 
and its components Bx,By,Bz transform as Bz ^ A2, 
{—By,Bx) ~ El under Ce^. They can couple only to 
valley-diagonal matrices. Indeed, for uniform electric 
and magnetic fields the effective hamiltonian must be 
invariant under translations, represented by the matri- 
. This requirement is not satisfied by any 
= Q,x,y, z. In partic- 



Ces e='=(27ri/3)A, 

matrix of the form EiA,., S,A„ 



ular, the term proportional to Bx^x^x + Byt^x^y, con- 
sidered in Ref. 45, must be absent (see Appendix B for a 
microscopic calculation) . 

These considerations enables us to write the hamilto- 
nian to the first order in the fields as 

^c-cm = -dz£z^ + [IxyiBx-^y ~ By^x) ~ i^zBz^z- (29) 

The first term in this hamiltonian represents the coupling 
to the z-component of electric dipole moment of tt elec- 
trons in each unit cell. The second term represents the 
coupling of the magnetic field to the in-plane component 
of the magnetic moment. These terms are forbidden for 
a suspended graphene sheet due to the symmetry with 
respect to reflection in the crystal plane, but may be al- 
lowed if this symmetry is broken due to the presence of a 
substrate. On the contrary, the third term, correspond- 
ing to the z-component of the magnetic moment of the 
unit cell, is allowed for a suspended graphene sheet as 
well (the z-component of the magnetic moment does not 
change sign under reflection in the crystal plane). 

Terms quadratic in £ and B correspond to polarizabil- 
ities of the unit cell. The corresponding hamiltonian can 
be written following the same lines as above. This is, 
however, beyond the scope of our interest. 



IV. RAMAN SCATTERING: GENERAL 
EXPRESSIONS 

In this section we derive the general expressions for the 
Raman scattering probability using the standard pertur- 
bation theory. 



A. Green's functions 

The second-quantized version of the Dirac hamilto- 
nian (10) reads as 



Hi = / d^rV'T(r) {-ivT. ■ V) V^r) 



(30) 



Since all energies we are interested in (~ 1 eV) are much 
higher than temperature, we set the latter equal to zero. 
The zero-temperature electronic Green's function, corre- 
sponding to hamiltonian (30), is given by 



G(p,e) = -« / hij{v,t)ij\Q,Q)\e'''P''+"-^(fYdt 

_ e -H vp • S 
e^ — {vp — io)"^ 



(31) 



where io is the inflnitesimal imaginary shift of the pole. 
T is the sign of the chronological ordering, and the aver- 
age is taken over the ground state of the system. 

Upon quantization of the phonon held based on hamil- 
tonian (20), the normal mode displacement operator and 
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the bare phonon hamiltonian become 



V- Kt^e"""" 



q 



5t g-i'i'^ 



^2NML0^{q) 



H, 



ph 



E 

q./j 



/ .^q.p 



"q,M"q,Ai 



E 



LxJ-ly 



(2^ 



(32a) 
(32b) 
(32c) 



The phonon Green's function is defined as 



^ 2NMu:M) ^Ta,(r,^)u,(0,0))e-'^-+-*d^rrft 



t^^ - [^A'(q) -«o]2' 



(33) 



This definition implies that each electron-phonon vertex 
corresponding to the second-quantized version of the in- 
teraction hamiltonian (23) contains, besides the coupling 
constant F^^ and the corresponding electronic matrix 
AjEj = (AE)^, a factor y/LxLy/[2NMLUf,{ci)]. Thus, the 
overall factor appearing in each vertex at small q is just 
Vy^X^/2, where A^ is the dimensionless coupling constant 
defined in Eq. (24). The factor L^Ly/N = V27aV4 is 
the area per one carbon atom. 

Upon quantization of the electromagnetic field the op- 
erator of the vector potential is expressed in terms of 
creation and annihilation operators oL ^, aq^^ of three- 
dimensional photons in the quantization volume V = 
LxLyLz with the wave vector Q and two transverse po- 
larizations f = 1, 2 with unit vectors eg ^i 



A(r) = 



Q,£ 



27rc 
VQ 



GQ/aQ^ie 



iQr 






(34) 

The photon propagator is defined analogously to the 
phonon one: 



. Q 



2ttc 



f (TAe{r,t)Ai{0,0)\e-'^''+'^'d^rdt 



2cQ 



n^ 



(35) 



so that each electron-photon vertex corresponding to 
the second-quantized version of the interaction hamilto- 
nian (28) contains the factor e{v/c), playing the role of 
the coupling constant, the electronic matrix (eQ,£ • S), 
and a factor y/2TTc/Q. 

Let us diagonalize each of the two electronic Green's 
functions, entering and leaving the electron-photon ver- 
tex, by transformation (11). Then, neglecting the change 
of the electronic momentum p upon emission of the pho- 
ton, we transform the electronic matrix in the vertex as 



„iS^i^pi/2piS„7r/4 

(P • eq. 



e""""^p''"e""«"' ' {eQj ■ S) e ""^fp/-e 



iE^V'p/2 -iS„7r/4 



[P X eq^^J 



^y, 



(36) 



where (pp is the polar angle of the wave vector p. Since 
we are interested in the interband transition, we need 
the T.y part of this expression, which tells us that the 
transition dipole moment is perpendicular to the electron 
momentum. 

The graphical representation of the Green's functions 
and vertices, introduced in this subsection, is shown on 
Fig. 7. 



B. Raman scattering probability 



Formally, n-phonon Raman scattering is a quantum- 
mechanical transition from the initial state with the crys- 
tal in the ground state and one incoming photon with the 
wave vector Qi„, frequency ujin — c|Qi„|, and polariza- 
tion Gin, into the final state with one outgoing photon 
with the wave vector Qowt, frequency LOout ~ c|Qo„i|, po- 
larization eout, and n phonons corresponding to normal 
modes /xi , . . . , ^„ with wave vectors qi , . . . , q„ . Denot- 
ing the ground state of the system by |vac), we represent 



vac) 



these two states as aj„|vac) and aLt^J^iqi • • • ^m„^„ , , 
Let us introduce the 5'-matrix S'(oo) and define the am- 
plitude of the n-phonon Raman scattering as 



A!^::Zi^;W^'S) 



T aout{t) bfj^^diih) 



^Mi.q- 



„(t„)5(oo)aL(0) 



jnt 



S{oo] 



dtY[, 



,iuJiti 



dU 



(37) 



Here the operators are taken in the interaction representation, and T represents the chronological ordering. 

The diagrammatic representation of the amplitude (37) in the leading order in electron-photon and electron-phonon 
coupling is shown in Figs. 8 and 9 for n = 1 and n = 2, respectively. The incoming photon line corresponds to pairing 
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of the operator aj^, the outgoing hnes - to pairings of the operators ciout, btii,qi- The loop represents the intermediate 
states of the electron-hole pair, summation over the momentum circulating in the loop and integration over energy 
should be performed. 

Let us separate the Green's functions of the scattering particles [only the positive-frequency parts of Green's 
functions enter, as shown by the {+) superscripts], the momentum-conserving i5-function, and explicitly introduce the 
permutations V of the phonon indices: 



AZ::Z{n;{c.,}iZi) = -^t'^M q„„o + ^^J -^Ti+)(Q„„„f7)n 



(27r)2<5( f^q, + Q„„, - Q,„ ) ^(-z)M;{^;}(r!;7'K}) 



Dl+\ci,,uj,)x 



(38) 



Then -~iA4 is given by the sum of all topologically inequivalent diagrams with amputated external lines. Equivalently, 
in the leading order '^'pi—i)Mv is given by the sum of all (n -I- 1)! connected pairings of electronic ■(/^-operators while 
electron-photon and electron-phonon vertices are held fixed. 
The transition probability per unit time is given by 

2 



lim — 

t^+oo t 



1 



where we used the relation 



/ n 



int d^ TT 



-iuit 



diOj 



27r J--'- 27r 



H 2 



(27r)2j[ ^qi + Qo„t-Qi. 



^m^|^;|K„,;7'KJ) 



lim — 

t^+00 t 



2it5{z) 



,(39) 



(40) 



The absolute dimensionless probability for the incoming photon to scatter with emission of n phonons of any kind 
is obtained by summing over all final states (here one should remember that a permutation of phonon arguments 
represents the same state) and multiplying by the photon attempt period Lz/c (at this point we also recall about 
the electron spin and multiply the matrix element by a factor of 2 which appears after tracing the fermion loop with 
respect to spin indices): 



/n = 



1 Lz 



11 / " 

■ ^ (L L Y ^ ^ 27r(5^a>^.(qO+c|Qo„t|-c|Q, 




H 2 



Ct^ 



2M 






(41) 



The square of the momentum (5-function is taken care of by the relation (27r)^(5(q = 0) = L^Ly, consistent with 
Eq. (32c). We can also pass to the spectrally resolved probability by inserting 1 = J dujout ^{'-^out — c|Qo„i|): 



dl„ 

dujoi 



^E 



-j^^S{c\Qout\ 



' ^out) 






) + ^out ~ ^v 



2M 



{q.} 



(42) 
At first glance, the matrix element A4)^'l docs not seem to depend on the outgoing photon wave vector QoMt, since 



'{q.} 



A^^^} 



the latter is negligible in comparison with electron and phonon momenta contributing to ■M)„ \, so the integration 

over the photon wave vectors gives just a>Q„j/(27r^c'^). However, Al| '! depends on the orientation of the polarization 
vector Coui, which, in turn, depends on the direction of Qout- Thus, we should consider the differential probability 
of emission into the elementary solid angle doout = sin 8 dQ d^, where the spherical angles & S [0, tt] and $ G [0, 2tt] 
parametrize the direction of Qowt- For this differential probability we can write 



47r 



din 



271" ^out 



dOn 



c 2'k'^(? (L^Ly 






{Mi}qi+---+qn 



2M 



{q.} 



^n (®m ' ^injl^out ' ^out) "r l-'n -'n ) 



(G'jri • ^out){^in ' ^out) < \^in ' ^out)\^in ' ^out) 



(43) 
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We stress that only the in-plane components of the polarization vectors participate in these scalar products. In writing 
Eq. (43) we have assumed the crystal itself to be isotropic, which is true as long as the calculation is done for the Dirac 
spectrum. As soon as the trigonal warping term in Eq. (12) is taken into account, the orienation of the polarization 
vectors with respect to the crystal directions will enter. The corrections due to the trigonal warping will be analyzed 
in Sec. VI B, and will be shown to be small. 

For each direction of Qout we can choose the basis of s- and p-polarizations: 



Gs = (— sin<I>, cos$, 0), Gp = (— cos 9 cos $, — cosGsin'I', sin 6). 



(44) 



Suppose the light coming out from the sample is collected by a lens within a cone with the aperture 2Qdet- Upon 
passing through the lens the polarization vectors change into 



(-sin<i>, cos<I',0), Gp -^ (-cos<I), -sin$,0). 



(45) 



so a linearly polarized detector oriented at an angle ^det will detect only those photons whose polarization before the 
lens was 



Gil = Gs sin(<i> - $ 



det) 



GpC0s(<I) - $ 



Averaging over the directions gives 



/„ = 



sin ede f d^ 



27r 





Odet 



det) 



(46) 



T^(o. . P* \J)- 4- (T\\ - T^\ in\^in! ^ \^tn) '^i- 



sin 6 dQ 
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-'n (,6m ' ^in) "ij ' \^n -'n I 



eL(eL)* + (eL)*e' 



(/ll + I^){1 - cosQdetf + 2/,|(8 - (1 + cosedet)') (e^ • e*„) + 



1 -" T^MQ n , „^..Q. ^3^l^- - m2 



+ ^«-^„)(8-(l+C0Se,et)') 



• Gdet ) 



(47) 



r 



where Gdet — (cos^^etiSin^detjO). Unpolarized detec- 
tion corresponds to adding the contributions from two 
mutually perpendicular polarizations Bdet and results in 



4 - 3 cos Qdet - COS'^ Qde 

12 



(/1+/^)(G„,.G*J. (48) 



C. Inelastic broadening 

As discussed in Sec. I A, when the number of emit- 
ted phonons is even, energy and momentum conservation 
can be satisfied in all elementary processes, represented 
by vertices on the diagrams for the Raman amplitude. 
As a consequence, the energy denominators of all elec- 
tronic Green's functions forming the electron-hole loop 
can be nullified simultaneously, and the integral over the 
internal momentum and energy diverges. To cure this 
divergence it is essential to include broadening of the 
electronic states. In other words, the infinitesimal imag- 
inary part —io in the denominator of the Green's func- 
tion (31) should be replaced by the actual broadening 
— Z7p — zImS(p, up), where S(p,e) is the electronic self- 
energy (the effect of ReS will be studied in Sec. VIII). 



ImS(p,fp) corresponds to emission of some excita- 
tions by the electron (hole) . One obvious candidate is the 
phonon itself, the corresponding contribution to S(p, e) is 
represented by the first term in Fig. 10. Besides phonons, 
an electron can emit other kinds of excitations, the most 
important contribution can be expected to come from 
emission of electron-hole pairs. Emission of electron- 
hole pairs in the undoped graphene must be impurity- 
assisted, ^^ while in the doped case electrons can collide 
without impurities involved, so the electron-electron col- 
lision rate strongly depends on doping. The propagator 
of electron-hole pairs is represented in Fig. 10 by the 
wiggly line. We do not need the explicit form of the 
propagator, being interested only in the contribution to 
IniE, introduced phenomenologically in the denomina- 
tor of G(p, e). The only property which is crucial for our 
consideration is that the spectrum of electron-hole pairs 
at each fixed wave vector is very broad as compared to 
Raman peak widths, as well as to 7p. 



The phonon contribution to E(p, e) (the first term in 
Fig. 10) can be calculated explicitly. Let us do it for the 
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case of E2 phonons: 
-*Sr(p,e) = 



F^ 



de' £p' V27^Dr{p~p',e-e') 

2^ (27r)2 4 2Afwr(p-p') 



X [S,A, G(p', e') S,A, + S^A, G(p', e') S^A,] (49) 

The two terms in the square brackets correspond to two 
phonon polarizations. Integrating over e' and neglecting 
the phonon dispersion, we obtain: 



S^(p,e) 



F^ 



27a2 



Uo. 



C'dc' 



Mturv^ 4 27r 7 £2 _ (-^/ + ^jp - io 



In 



An 



ujr 



-{e^^e) 



(50) 



where ^max is an ultraviolet cutoff of the order of 
the electronic bandwidth. If we consider scattering on 
E[ phonons (Ai representation of C^y in the K point), 
we obtain the same expression as Eq. (50), but with Fr 
and ujr replaced by Fk and lok- The dimensionless cou- 
pling constants Ar and Xk are defined in Eq. (24). 

To conclude this section, we note that Raman scat- 
tering rate can be viewed as the imaginary part of the 
incoming photon self-energy; indeed, any diagram for the 
Raman matrix element, e. g., the first one in Fig. 9, can 
be obtained by cutting a photon self-energy diagram in 
two, as shown in Fig. 11. Dressing the electronic Green's 
function with self-energy corrections corresponds to go- 
ing to higher orders of perturbation theory and taking 
into account diagrams of the type shown in Fig. 12. One 
may ask, what other diagrams may appear in the higher 
orders of perturbation theory, and how important they 
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FIG. 7: Graphical representation of electron, phonon, and 
photon Green's functions and vertices. 





(a) 



FIG. 8: Diagrams for the one-phonon Raman amplitude in 
the leading order. 



• Diagrams where the two electronic loops are con- 
nected with more phonon lines (Fig. 13) describe 
Raman scattering on the corresponding number 
of phonons. If the loops are connected with an 
electron-hole propagator, this corresponds to the 
contribution of electron-hole pairs to the Raman 
spectrum. As discussed in Sec. I A, the resulting 
spectrum is broad and featureless, so we are not 
interested in these processes. 

• One can also insert phonon and electron-hole prop- 
agators as shown in Fig. 14. Such diagrams repre- 
sent vertex corrections, their effect is analyzed in 
Sec. VIII. 

• All diagrams involving more electron-photon ver- 
tices have an extra smallness in the parameter 
e^/c « 1/137, and are neglected. 

• Inserting electronic loops into phonon propaga- 
tors describes the phonon frequency shift due to 



the continuum of electron-hole pairs. They have 
been studied before, both theoretically^^'^^"'** and 
experimentally.^' ^°''^°'^^ However, shift and broad- 
ening of the phonon states are not important for our 
calculation, as we are interested in the frequency- 
integrated intensities of the Raman peaks. They 
are determined by the total phonon spectral weight, 
which is not changed by shift and broadening. 



V. ONE-PHONON RAMAN SCATTERING'"' 

Since the photon wave vector is negligibly small com- 
pared to all other scales in the problem, the in-plane 
momentum conservation requires the emitted phonon to 
belong to the F point. This fact explains the small width 
of the peak at 1580 cm~^. The matrix element for the 
one-phonon process is given by two diagrams in Fig. 8: 
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(b) 




FIG. 9: Diagrams for the two-phonon Raman amplitude in 
the leading order. 



-zS(p,e) = 



+ 




FIG. 10: Electron self-energy due to the interactions with 
phonons and with electron-hole pairs. The propagator of 
electron-hole pairs is shown by the wiggly line, its explicit 
form is not important for our calculation. 




FIG. 11: Raman scattering rate as the imaginary part of the 
photon self-energy: diagram (a) on Fig. 9 for the Raman ma- 
trix element is obtained by cutting the photon self-energy 
diagram in two, as shown by the vertical short-dashed line. 




FIG. 12: A diagram corresponding to self-energy insertions 
in the electronic Green's functions. 
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FIG. 13: A diagram describing emission of extra excitations. 




FIG. 14: A diagram representing vertex corrections. 
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27re2i;2 



iy/LOinUJout y 2NMujt{0) 



J^xJ^yi'Y 



de d^p 



Tr {(e„S) G(p, £_ - ur) T,,y G(p, £_) (e:„,S) G(p, e^ 



2^ (27r)2 rx4 ^^^■'""' ^ ^^' ^'- "^ ' ^"'^ 

+ G(p, e_) (e:„,I]) G(p, e+) T,,, G(p, e+ + c^r) (e^S)} , (51) 
I 



where we denoted T^ = A^Ey, Ty 
ajo„t/2 and wr = wr(q = 0). 



-A,S, 



e± 



e± 



At first glance, elementary power counting (we set 
e ^ vp and note that each Green's function G ~ 1/e) tells 
us that the integral is logarithmically divergent at high 
energies. However, let us recall the continuous symme- 
try (15) of the Dirac hamiltonian (10). A rotation (15b) 
of the dummy integration variable p by an arbitrary fixed 
angle tp cannot change the value of the integral. Thus, 
the integral must not change if each Green's function 
is replaced by G ^ g-iT,,ip/2 Q^i^,ip/2^ Application of 
these e^^^^v/2 to the vertices is equivalent to the inverse 
rotation (15b) of the polarization vectors ei„,e*^^ and 
of the matrices E^; , Ey . This gives us a cubic combina- 
tion of smip,cosip which has no (/s-independent terms. 
This may be viewed as conservation of the z-component 
of the angular momentum: indeed, transformation prop- 
erties under rotations Coov for both the photon polar- 
ization vectors ei„ , Bout and the E2 phonon displace- 
ments UE2x,UE2y correspond to the angular momentum 
m — ±1. As a result, for the Dirac spectrum the one- 
phonon matrix element (51) must vanish. 



In fact, Eq. (51) gives zero even prior to the p- 
integration. We just notice that 

-zA^G(p,e)zA, = G(p,e), (52a) 

-zAa^SiA^ = S, (52b) 

-iA,(A,S)zA, = -A,S. (52c) 

The matrix —iA^ is the combination of the matrix A^jSz 
(the full G2 rotation) and e~'^^''"/^ (rotation by tt within 
each valley). This symmetry relates the spectra at the 
points K -I- p and K' -f p. 

As a consequence, in order to describe the one-phonon 
Raman peak at 1580 cm"-*^, one has to go beyond the 
leading order in the small-p expansion of the hamilto- 
nian [i. e., to take into account H2 from Eq. (12)] and 
interaction vertices. Since Aa;-symmetry (52a) is broken 
by the second term of the hamiltonian H2{p) [Eq. (12)], 
the result will be different from zero already at the next 
order in p. Then an additional power of p appears in 
the integrand of Eq. (51), so the divergence at the up- 
per limit becomes linear rather than logarithmic. This 
means that the small-p expansion is inapplicable (i. e., all 
of its terms give the contribution of the same order) , and 
the whole Brillouin zone is responsible for the 1580 cm~^ 
peak. We stress that this statement is valid in the clean 



22 



limit, when impurity scattering can be neglected. Impu- 
rity scattering can allow the one-phonon process in the 
leading order and make the integral convergent in the 
upper limit. ^ In the clean limit the proper tool for the 
calculation is thus not the effective low-energy electronic 
theory, but ab initio band structure methods, which are, 
of course, beyond the scope of the present work. For 
low photon energies we are interested in, we can simply 
add a term to the system hamiltonian, corresponding to 
a direct photon-phonon interaction. The form of this 
hamiltonian is fixed by the Cq^ symmetry (Ei x Ei x E2 
contains only one Ai representation): 



H, 



cm— ph 



VR 



/ d^^ [i^x - £yhE2V - '2£x£yUE2x] , 



(53) 

where £x,£y are the in-plane cartesian components of the 
electric field vector £ = —{l/c)dA/dt, e is the electron 
charge, and ur is the unknown constant of the dimension- 
ality of velocity. Since it originates from the electronic 
TT-band, its magnitude is roughly given by the product of 
the electronic bandwidth and the lattice constant, i. e., it 
should be of the same order as the electronic velocity v. 
This hamiltonian leads to the following expression for 
the one-phonon Raman matrix clement: 
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271 y/LUinUJout 
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(54) 



-■in,x'^out,x 



The resulting intensity does not depend on the polariza- 
tion of the detector. The sum over the two polarizations, 
calculated according to the prescription of Sec. IV B, is 
given by 

/r = Stt 



w2c2 4Mwr(0) 
4 — 3 cos Qdet — cos'^ 6 



'Jdet 



'Jdet 



(er 



(55) 



Furthermore, Eqs. (52) lead to vanishing of the Ra- 
man amplitude for any odd number of T phonons. This 
property resembles Furry's theorem in the spinor quan- 
tum electrodynamics. ^° There is, however, a difference: 
Furry's theorem holds for the sum of two diagrams con- 
taining an odd number of photon lines and differing by 
the direction of the electronic loop which cancel each 
other, while in Eq. (51) each of the two terms vanishes 
separately. We also note that K phonons can be emitted 
only in pairs: the excited electron-hole pair should switch 
valleys an even number of times in order to annihilated^ 
The difference from the one-phonon process is that the 
integral in the matrix element converges at large k even 
after the next term in the small-A: expansion have been 
picked up. Thus, the three-phonon process can be de- 
scribed within the low energy theory, but the corre- 
sponding amplitude will contain an additional smallness 



to. 



la/c. 



VI. TWO-PHONON RAMAN SCATTERING 

A. Calculation for the Dirac spectrum 

First, let us focus on the 2K peak at 2700 cm^^, which 
corresponds to emission of two scalar phonons (Ai in 
terms of the C^v symmetry, or E[ in terms of the Cg„ 
symmetry) from the vicinity of the K and K' points. 
The integrated intensity of the 2K peak is given by 



47r 
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2K 



dOn 



2^ <ut 1 
c2 27r2c2 2 



(27r) 



\2M 



K\ 



(56) 



The factor 1/2 in front of the sum eliminates double 
counting in the summation over the final states (phonon 
permutations), the factor of 2 comes from the sum over 
Ai and Bi modes (no cross-terms arise, as they would 
yield a traceless combination K^-l^y), the factor of 2 in- 
side the square takes care of the spin degeneracy. 

The diagrams for the two-phonon matrix element are 
shown in Fig. 9. Only the first diagram corresponds to 
the fully resonant process, described in Sec. I A, the other 
two give a contribution, smaller by a factor ^ /ujk- Thus, 
the matrix element is given by 



M 



K 



27re2v2 L^LyF^ 



de d^p 



y/i^tn^out 2NMcjK{q) J 2tt (27r)2 

X Tr^ {(e.„ • S) G(p, e_) A,S,G(p + q, e'_) x 

X (e:„, • S) G(p + q, eV) A,S,G(p, e+)} + 
+ (q^-q), (57) 

where e± = e ± ujin/2, e'_j_ = e ± ujout/'^- We diagonalize 
the Green's functions by the unitary transformation (11) 
and neglect the off-resonant contribution: 



G(P, e,) = L/t f±±^ L/p « C/t (i±|^ C/p , 

(58) 
where Cp = v\p\ - ijp, Up = e'^y^/^e'^'"^^/^, and tpp = 
arctan(pj,/pa;) is the polar angle of the vector p. The 
unitary matrices rotate the vertices as 

Up Sj; U^, = T.^ cos ^ - ^y sm ^ (59a) 

Up E, Ul, = S, sin ^^±^ + S, cos ^P±^(59b) 

C/pE,C/t, = zlsin ^P'^P' - S.cos ^P~^P' ,(59c) 

Up (e„i • E) C/jl, = (ep • ej„) S^ + [ep x e,n]z Sy,(59d) 

where Gp = p/|p|- We are interested in the Sj, term in 
Eq. (59d), which corresponds to interband transitions, 
and also tells us that the transition dipolc moment is 
perpendicular to the electron momentum. In Eq. (59c), 
the rotated phonon vertex, we need the intraband term 
ex 1. Evaluation of the trace gives 

2 [ep X ei„]^[ep+q x e^^J^ sin^(y)p+q/2 - ipp/2) 

(6_+ep)(£'-+ep + q)(e;-ep+q)(e+-ep) ■ ^ ' 
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We can rewrite each of the factors in the denominator of spend in the corresponding intermediate states by the 

electron and the hole, respectively. 

Let us denote pa = uji„/{2v), pi = Wo„t/(2w), q = 
Q — Po ^ Pi, and let P||,p_L be the components of the 



this expression as 
1 
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(61) 



Integration over de/{2'!r) produces 6(to + ti — to — ti). 
Thus the times io, ^i and to, ii can be interpreted as times 



J 



deviation p = p + po^q along q and perpendicular to q, 
respectively. We expect the deviations to be small, so we 
approximate 
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(62a) 



(62b) 



and Gp w -Gq, Gp+q w Gq, (/3p+q/2 - (/)p/2 K. 7r/2. Inte- 
gration over momentum deviation p is performed as (we 
denote 7po + 7po+q = 27 for brevity): 



dto dti dto dti 5{to 



, ti - io - ii) e~™'^(*i+*i' / — e™J'ii(*o+*o"*i~*i) 



2tt 



X / !^(,-bpo+'^Pl/C^Pomto+io)-hpa+^+ivpl/{2pi)]{U+h) ^ 1 ^ / CJ,„tJo~ 1 
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(63) 



r 



which gives the matrix element: 



M^ = 



2ne^v^ L^LyF^ 1 

y/ujin + ujout 2NMujK{q) Sv"^ 

Z[Gq X Gi„J^[Gq X G^ 



{vq - 2i7)3/2 



+ (q^-q). (64) 



This expression describes a strongly peaked q depen- 
dence of the matrix element Al(q) around the value 
qbs = {tOin + <^out)l(2v), corresponding to backscatter- 
ing of the electron and the hole by the phonons. The 
width of the peak is determined by the electron life- 
time: Sq ~ j/v. Such sharply peaked dependence cannot 
be derived from pure symmetry considerations, which 
just prescribe vanishing of the phonon matrix element 
for the forward scattering.^ It is a consequence of the 
fully resonant character of the two-phonon Raman scat- 
tering, which should be contrasted to the double-resonant 
impurity-assisted single-phonon scattering. For the lat- 
ter case, the scale determining the width of the peak is 
6q ^ ujin/v [see Eq. (2) of Ref. 51]. 



The origin of such peaked dependence lies in the quasi- 
classical nature of the electron and hole motion. Namely, 
the incoming photon creates and electron with momen- 
tum Po, moving with the velocity vq = 9Re^p/9p|p=pj,, 
and a hole with momentum — Po, moving with the ve- 
locity — Vq. After a time to the electron emits a phonon 
of momentum — q, and after a time to the hole emits a 
phonon of momentum q. Afterwards, at time to + '^i = 
to + ii they recombine and emit a photon. In order to 
recombine, they must meet at the same spatial point 
(as prescribed by the spatial (5-function resulting from 
the integration over the deviations p) with opposite mo- 
menta pi = Po + q and — pi (as prescribed by the mo- 
mentum conservation) . It is possible only if the velocity 
vi = 9Re^p/(?p|p=pj^ is directed oppositely to the ini- 
tial velocity vo. The deviation of the typical scattering 
angle (p from tt is ji/? — 7r| ^ \/l I'^in ^ 1, determined by 
the quantum diffraction [note that we had to expand the 
energy to the second order in p^ in Eqs. (62)]. 

Proceeding with the calculation, we substitute Eq. (64) 
into Eq. (56). Angular averaging gives 
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-TT^ |[eq X e^n]z[eci x e*„t]^|^ = - [(ei„ • eo„t)(e*„ • e*„j) + (ej„ • e*„t)(e*„ • eo„t) + (ej„ • e*„)(eo„t • e*„t)] , (65) 

ZTT 8 



SO that /2jf = 3/^, as seen from Eq. (43). We stress that only the in-plane components of the polarization vectors 
participate in these scalar products. Evaluating the q integral and using Eq. (47), we obtain the final result: 
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(66) 



where qi,s = {(^in + '^o«t)/(2w) is the phonon wave vector 
corresponding to the backscattering, and V27a^/4 is the 
area per carbon atom. Eq. (8) of Ref. 28 corresponds'"'^ 
to normal incidence (|ei„| = 1), collection in the full solid 
angle 47r, and summation over the two orthogonal direc- 
tions of e^et, which makes the second line of Eq. (66) 
equal to 1/3. 

Eq. (66) shows how the dominant role of backscat- 
tering manifests itself in the polarization memory. In- 
deed, linearly polarized light preferentially excites elec- 
trons and holes with momenta perpendicular to the elec- 
tric field vector. After the phonon emission these mo- 
menta change to the opposite, and the photon emitted 
after the annihilation has a preferred direction for the 
polarization, perpendicular to the electron and hole mo- 
menta. Quantitatively, Eq. (66) gives the ratio of inten- 
sities for edet\\ein and e^et -L ei„ to be 3 at Odet -^ and 
23/9 at Qdet = 7r/2 (collection into the solid angle 2tt). 

Let us now turn to the 2T phonon peak at 
3250 cm^^, corresponding to emission of two pseudovec- 
tor E2 phonons from the vicinity of the F point. Its 
integrated intensity can be calculated analogously: 
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(67) 



As in Eq. (56), the factor 1/2 in front of the sum elim- 
inates double counting in the summation over the final 
states (phonon permutations), the factor of 2 inside the 



square takes care of the spin degeneracy. Here, however, 
the summation over the two phonon modes is less simple. 
First, let us neglect the phonon dispersion, so that the 
longitudinal and transverse phonons are degenerate. The 
matrix element is given by 



M'^ 



27re2w2 



^x^y^r 



de d^p 



y/nJ~IJ^t 2NMLJT J 27r (27r)2 

X Tr^ {(e„, • S) G(p, e_) T, G(p + q, e'_) x 

X (e:„,.S)G(p + q,6V)r,G(p,e+)} + 

+ (q^-q), ^ (68) 

= A^Ey, Ty — — A^Sa;. Thc unitary transfor- 
mation (11) rotates the isospin as given by Eqs. (59a), 
(59b). Thus, the trace will be given by an expression, 
analogous to Eq. (60), but with the replacement 



where T-r 



sin2((^p+q/2 - (/?p/2) -^ -e'e-' 



cos 



y^p + yp+q 



All arguments leading to the dominance of backscattering 
remain valid, so we obtain e — — eq which means that 
only the longitudinal phonons are emitted. Summation 
of the probability over the two phonon polarizations gives 
(e^)4 + {eyf + 2(e^)2(e«)2 = 1, so the intensity of the 
2F peak is given by one half of the expression (66) with 
replacements Fk -^ Fy, ojk -^ ^r,L- 
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2\cJ c2 872 \Mv''Lur,L{qbs) 4 
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I — 



8- {l+COsQdetf .. ^|2 
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(69) 



B. Effect of trigonal warping and electron-hole 3250 cm ^, respectively. According to the data of Ref. 6, 

asymmetry 



In the previous subsection we have calculated the inte- 
grated intensities of 2K and 2F peaks at 2700 cm^^ and 
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hn/hr ~ 20. At the same time, for the corresponding 
electron-phonon coupling constants the nearest-neighbor 
bond-stretching approximation gives Fk/Fy — 1, which 
agrees with DFT calculations of Ref. 29 with the preci- 
sion of 1%. Then, what is the origin of such huge differ- 
ence in the intensities of the two peaks? In this section we 
consider the effect of electron-hole asymmetry and trig- 
onal band warping on the intensities of the two-phonon 
Raman peaks with the purpose to check whether the trig- 
onal warping can explain the observed large difference of 
intensities I2K and /2r- 

Typically, one neglects corrections to the Dirac spec- 
trum u|p|, arising from the quadratic term i?2(p) in the 
electronic hamiltonian, given by Eq. (12), as they are 
smaller than u|p| by a factor pa <^ 1. However, accord- 
ing to the results of the previous subsection, two-phonon 
scattering is sensitive to the directions of electronic veloc- 
ities and momenta on the angular scale Sip ~ y^j/ujin- 
This means that effects of (i) phase mismatch between 
the electron and the hole, introduced by the first term 
in Eq. (12), and (ii) non-coUinearity of velocity and mo- 
mentum, introduced by the second term in Eq. (12), may 
become important already when H2{p) ^ 7, which hap- 
pens at a smaller energy scale than H2 (p) ~ if 1 (p) . This 
means that, while including H2{p), we still can neglect 
higher-order terms of the expansion (-ffs, H4, . . .). 

Thus, we repeat the calculations of the previous sub- 
section, taking into account the H2(jp) term only in 



the denominator of expression (60) (the numerator is a 
smooth function of p, so corrections to it will be small as 
pa indeed). The quadratic term (12) in the hamiltonian 
modifies the electron dispersion as 
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dp 
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- = cosSipp, 



p.:, 



(70a) 
(70b) 

(70c) 
(70d) 



3P^ - ^plPy 



po 



-3sin3(/3p, (70e) 



where the sign of the ao terms is "-I-" in ^p (electron 
dispersion) and "— " in ^p (hole dispersion), and the sign 
of the a^ terms is "-f" and "— " for K and K' valleys, 
respectively. 

For a given direction Gq we choose po = — po^q, Pi ~ 
piGq, with Po and pi such that Re^pj, + Re^p„ — ujim 
Re^pj + Rel^pj = ujout, and denote q = |q| - Po - Pi, 
p = p — Po- The energies of the intermediate states to 
be substituted in Eqs. (61) at the corresponding times, 
can be approximated as 
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(Po •p) + 2a3Z^_L(Po)[Po X p]^ + 7— giPo X p]l, 
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(71a) 



(Pi • P) T 2a3i^_L(pi)[pi X p]^ + -—[pi X p]j + vq, (71b) 



(po • p) + 2a3U^{pa)[po x p];^ + 7r^[po x p]^. 



(Pi • P) T 2a3Uj_{pi)[pi X p]^ 



(71c) 



V 

2Po 

^[pixp]l+vq. (71d) 
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As before, the upper (lower) sign corresponds to emission of scalar E[ (pseudovector E2) phonons, accompanied by 
intcrvalley (intravalley) electron scattering. Then instead of Eq. (63) we have 
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(72) 



Let us denote 7pQ + 7pj — 27, and introduce T — Iq + Iq — ti + ti (the latter equality follows from the longitudinal 
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spatial (5- function in the leading order), r = to — to = — ti + ti, and rewrite the integral as 
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where we also denoted 

Aeh = ao{pl -pf). 
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Squaring the matrix element, and performing the final 
integration, we obtain Eq. (66) with the replacement 
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eh 



In 



7^ + AL/4 
7^ ■ 



(75) 



The meaning of this replacement is that when the 
electron-hole asymmetry becomes greater than the level 
broadening, it starts to play the main role in restrict- 
ing the energy denominators from below. Numerically, 
ao(l eV)^/w^ ~ 0.1 eV (see, e. g., Ref. 19), so the rela- 
tive correction to Eq. (66) for small A^h can be estimated 

as -(1/2)A2^/(27)]V2 10-\iJ^n/2J)^. The total 

electronic broadening 27 was measured by time-resolved 
photoemission spectroscopy to be 20 meV in Ref. 20 and 
25 meV in Ref. 21 (all values taken for e = Win/2 — 1 eV). 
A recent ARPES measurement gives a significantly larger 
value for 27 ^^ 100 meV (Ref. 22), which agrees better 
with Eqs. (50), (25). 

The effect of the trigonal warping term (the only one 
sensitive to the difference between intervalley and in- 
travalley scattering) turns out to be small in the param- 
eter a^ujin/v'^ [in Eq. (73) warping enters through Aj^, 
which results only in a small shift of the integration vari- 



able q]. For uJin = 2 eV, a^ 



-va/A (tight-binding 



model) , we estimate the relative contribution of the warp- 
ing term as A±/uJin ~ 5 • 10""*, 5 • 10~^ for intervalley 
and intravalley scattering, respectively. Thus, trigonal 
warping cannot account for the observed ratio I2K/ l2r- 



VII. FOUR-PHONON RAMAN SCATTERING 

A. Resonant manifold 

We will calculate only the intensity of the peak which 
we will call AK^ the double of the 2K peak at 2700 cm~^. 



The AK peak corresponds to emission of four scalar 
phonons {Ai in terms of the C^^ symmetry, or _E( in 
terms of the C'^^ symmetry) from the vicinity of the K 
and K' points. We will not consider other four-phonon 
peaks, as their intensity is smaller. 




FIG. 15: The only fully resonant diagram for the four-phonon 
Raman amplitude in the leading order. 



The integrated intensity of the AK peak is given by 
dliK 27r uJl^^ 1 8 



47r- 



dOn 



c2 27r2c2 [L^Lyf 4! 



E 



\2M- 



qi:-q2,q2.q 



,q.,qj'.(76) 



qi-l-q2=qi+q2 



The factor 1/4! in front of the sum ensures the proper 
summation over the final states (phonon permutations), 
the factor of 8 comes from the sum over Ai and Bi modes 
(one can have all four phonons Ai, all four Bi, or 
4!/(2! • 2!) = 6 combinations of two Ai and two Bi 
phonons, all other combinations yielding traceless prod- 
ucts of Kx and Ay), the factor of 2 inside the square takes 
care of the spin degeneracy. The only fully resonant di- 
agram for the four-phonon matrix element is shown in 
Fig. 15. The matrix element is given by the sum of 
4! = 24 permutations of the emitted phonon wave vectors 



-qi, -q2, q2, qi [we denote ujk^o) = Wq 
compactness] : 



for the sake of 



27 



■'^-qi,-q2,q2,qi 



27re2w2 [L^LyFf^/{2NM)Y f de d^p 



\/i^in^out ^WqiWqjWqjW, 



q2 



27r (27r) 



X Tr^ { (ei„ • S) G(p, e - Wj„/2) E^G(p + qi , e - Wj„/2 + tjq^ ) S^G(p + qi + q2 , e - a)j„/2 + tjq^ + cjq^ ) x 

x(e*„( • I])G(p + qi +q2,e + Wi„/2- Wqj - WqJ S2G(p + qi, e + a)j„/2 - cjqj S2G(p, e + a)i„/2)} + 
+ (23 other permutations of q's). (77) 

Again, we switch to the basis of Dirac eigenstates according to Eqs. (58), (59), and evaluate the trace: 

[ep X e,„]^ sin(yp/2 - ipp+^j2) sin(yp+q,/2 - fp+q,+qj2) 



Tr{...} = 2 



(e - UJin/2 + ^p)(e - UJin/2 + UJ^^ + Cp+qi)(e - ^m/2 + Wqj + UJq2 + Cp+qi+q2) 

[ep+qi+q2 X Kut\z sin(yp+q,+qj2 - ipp+qj2) sin{ipp+qj2 - ipp/2) 

(e + tJi„/2 - Wqi -UJq,^ -^p+qi+qJ(e + W,„/2-tJq, - ^p+q J (c + Wi„/2 - ^p) 

I 



(78) 



Just hkc for the two-phonon scattering, the dominant 
contribution to the integrals will come from those mo- 
menta which make small all the factors in the denom- 
inators of the above expression. The real parts of the 
electron and hole dispersion, Re^p, Re^p, are given by 
Eqs. (70a), (70b), and Im^p ::= Imfp = ~jp. 

Four phonon wave vectors — qi, — q2, q2, qi, such that 
qi + q2 = q2 + qi , will be said to satisfy resonance condi- 
tions if exists a vector po such that the following equali- 
ties hold: 

J^S^po — „ ) 



RS^po+qi 
R^spo+qi 
Re^poH 

R6 4po+qi+q2 



Wq 



^qi+q2 



'"qn 






- UJc,, - LOr, 



^qi 



q2- 



(79a) 
(79b) 
(79c) 
(79d) 
(79e) 



The structure of the manifold defined by these condi- 
tions essentially depends on (i) whether we take into ac- 
count the electron- hole asymmetry, or just set Re^p — 
Re^p = ^'IpI, and (ii) whether we take into account the 
phonon dispersion, or just set Wq = wq- Indeed, subtract- 
ing Eq. (79e) from Eq. (79d), we obtain 



'"qi 



^q2 



ReCpo- 



-qi+q2 



^qi + ^q2 



RcCpo^ 



>po+qi+q2i 
(80) 
which in the phonon-dispersionless case represents either 
an identity if ^p = ^p, or can never be satisfied if ^p ^ ^p, 
thus leading to a finite energy mismatch; only in the 
phonon-dispersive case it represents a non-trivial equa- 
tion. As will be seen in Sec. VII F, the quantitative con- 
dition for the phonon dispersion to matter is that the 
phonon group velocity Vph > f max{7, Ae^}/aji„, where 
7 is the typical value of the electron broadening, and 
Ag/i is the typical value of the electron-hole asymmetry, 
defined analogously to Eq. (74b). 

First, let us focus on the electron-hole symmetric and 
dispersionless case, Re^p = Re^p = w|p|, cjq = ojq, when 




FIG. 16: Structure of the resonant manifold of electron and 
phonon momenta for four-phonon Raman scattering in the 
approximation Re^p = Re^p = u|p|, tjq = tJo- 



the structure of the resonant manifold (79) is the sim- 
plest. The resonant manifold can be parametrized by four 
polar angles (/?o, <^i, <^i, </'2, which determine the positions 
of the four momenta po, pi = po + qi, Pi = po -I- qi, 
and p2 = Po + qi + q2 = Po + qi + q2 on the three 
circles i;|po| == tjj„/2, w|pi| == w|pi| == tJm/2 - wq, 
i;|p2| — LOinl2 — 2ijjQ (Fig. 16). Since the total num- 
ber of independent phonon variables is 6, there are two 
more variables besides the four angles, which correspond 
to deviations of the phonon momenta from the resonant 
manifold. There is a freedom of choice for them, and we 
choose them to correspond to stretching of the vectors 
Pi,Pi (denoted by w,t(;, respectively): 

(5qi = -[e^ X po](5(^o + {^z x Pi]<5¥^i + Pi(5w;, (81a) 

^<\2 = \^z X P2](5(^2 - \Qz X Pl]<5¥'l - pi<5w, (81b) 

i5qi = -[e^ X po](593o + [e^ x pi](5(^i -I- pi(5w, (81c) 
^^2 = [e^ X P2]5(^2 - [ez X pi](5<^i - pi(5w, (81d) 
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The q integration measure is transformed as 

(fq_i<fq2d^q.i = IpiHpiP |[Po x pz]^! x 

X difio dip2 difii dipi dw dw. (82) 

Thus, 6 independent integration variables, parametrizing 
the final state, have been separated into two groups: four 
angles, determining the position of the phonon momenta 
on the resonant manifold, and two deviations w, w from 
the manifold. The dependence of the integrand on the 
first group of variables is smooth, and we will call them 
slow variables. The deviations from the manifold, on the 
contrary, suppress the integrand dramatically, so we will 
call them fast variables. 

Let us now consider other 23 permutations of the 
phonon wave vectors — qi, — q2, ^2, qi, which have to be 
added in order to obtain the correct value of the matrix 
element Ai. We need to check whether contributions 
to M from different permutations can be large at the 
same time (i. e., for the same values of the integration 
variables — qi, — q2, q2j Qi)- If this is not the case for all 
23 permutations, then the summation over all permuta- 
tions in Ai would just cancel the factor 1/4! in Eq. (76). 

First, we note that if three momenta Po, Po + Qi, Po + 
qi + q2 satisfy the resonance conditions (79) (i. e., they 
lie on the corresponding circles) for some po, then the 
momenta Po, Po + q2, Po + Qi + ^2 do not satisfy these 
conditions for any Pq. Indeed, the vector qi + q2 can 
be placed on the circles only in two ways, which leaves 
only two choices for the momentum pg: Pg = Po, or the 
symmetric one. For none of them does Pg + q2 lie on the 
circle, unless the configuration of momenta has special 
symmetries. This means that the corresponding contri- 



bution to the matrix element is small, so the permutation 
qi ^-> q2 should be discarded. By analogous argument we 
can discard all 14 permutations which leave at least one of 
the q's in place (this immediately leaves only two choices 
for Pq which can be inspected). 4 cyclic permutations 
are eliminated because the circles have different radii. 
The permutation qi <^ q2 , qi ^^ q2 is eliminated by the 
very first argument. Thus, we are left with four permuta- 
tions. One is -qi, -q2,q2,qi -^ qi,q2, -q2, -Qi, audit 
does satisfy the resonance conditions (79) if one chooses 
Pg = — Po- The other three are obtained from it by swap- 
ping the first two momenta, the last two, or both, and 
thus do not satisfy the resonance conditions. The per- 
mutation -qi,-q2,q2,qi —" qi,q2,-q2,-qi is nothing 
else but the result of reversal of the electronic line direc- 
tion in the loop on Fig. 15. As seen from Eq. (78), this 
permutation gives exactly the same contribution as the 
original one. Thus, their interference results in an addi- 
tional factor of 2, besides cancellation of 1/4! in Eq. (76). 



B. Quasiclassical representation 

Now let us pass to time representation for the matrix 
element M by rewriting each factor in the denominator of 
expression (78) analogously to Eqs. (61). We introduce 
three time variables for the electron and three for the 
hole, denoted by to, ii, ^2 and to, ?i, ?2, respectively. Next, 
we introduce the deviation p = p — Po from the point po 
fixed by the resonance conditions, and expand each factor 
in the denominator of expression (78) as 



J 



*" 1 I c 



~2~ 



- t^n, + i^n 



q2 I Sp+qi+q2 
Win 



Cp 



-^ ^qi ^ 4p+qi 



Cp+qi 



vop- «7o - aoPo) 

viP - Z7i -I- wvpi + Sujq^ - aoPi, 

V2P - J72 + Sujq^ + Sujq^ - Uapj, 

-vqp + ijo - aapl, 

-vip + i7i - wvpi - Sujq^ - aopj, 



+q2 



-V2P + ij2 - Su;^^ - (5wq 



aoP2, 



(83a) 
(83b) 
(83c) 
(83d) 
(83e) 
(83f) 



where v^,— v^ and 7i,7i are the velocities and the damping rates of the electron and the hole, respectively, in the 
ith intermediate state. In Eqs. (83) we have also taken into account the electron-hole asymmetry and the phonon 



dispersion 5ujc. 



OJa 



luq, both assumed to be weak: aoPg j^ 2 ^ ^o, I'^'-^ql ^ '-^o- Within this approximation we 



take Vg = Vg and V2 = V2- 

The numerator of Eq. (78) is a smooth function of momenta on the scale "f/v, so it can be taken at p = 0, i. e. on 
the resonant manifold. Then the integration over (i^p/(27r)^ gives a spatial (5-function: 



_fp 



-ip(voto+Viti+V2t2+Voto+Vltl+V2t2) 



<5(voto + Viti + V2t2 + Vgtg + Viti + V2t2) 



so in order to recombine, the electron and the hole should meet at the same spatial point. As a result, the denominator 
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of Eq. (78), integrated over de/{2n) and (i^p/(27r)^, is rewritten as 

e^** (5(to +ti+t2 — iQ — ti— h) S{voto + viti + V2t2 + vo^o + viti + V2t2) dto dii (it2 dio dii dh, (84a) 

2 

r] = i [wvpiti +wvpiti +dujq^(ti +^2) + '^t^q2*2 + <5'^qi(^i + h) +5ujqj2] +»ao^Pi(ij - ij) + 



'^{itU + jiU). 



i=0 



(84b) 



r 



C. Integration over deviations 

The part of the summation over final phonon states in 
Eq. (76), corresponding to the integration over the devi- 
ations WjW, can be performed exphcitly. Let us open the 
(5- functions in Eq. (84a) choosing tg, ti, fi as independent 
variables: 



to + to — 
t2 + h = 

t2,t2 — T [ to 



[V2 X (Vi^i +Vlii)]z 

[vo X V2]2 
[vq X (Vjtl +Vltl)]z 

[V2 X Vo]z 

ti - h 



[(vo TV2) X (viti +viti)]j 



2[V2 X Vo]z 

The Jacobian of this transformation is given by 

dto dt2 dt2 5{to + ti + t2 - to - h - h) x 

X (5(vo(to + to) + V2(t2 + h) + viti + viti) 
1 



(85a) 
(85b) 

(85c) 



2|[vo X V2]z\ ' 

The detuning phase from Eq. (84b) can be written as 



denote by O, besides the condition Iq > 0, is determined 
by the inequalities 



to < 
to < 



[V2 X (Vi^i +Vlti)]z 
[vo X V2]2 

h - ti [(v2 - Vq) X (viti + viti) 
2 2[vo X V2]^ 



(88a) 



ti-h , [(v2+vo)x (Viti+Vi^i)]^ 
to > 7; 1 TTr 1 .[OliC) 

2 2[vo X V2]^ 



Consider integral (84a). Squaring its modulus, and 
integrating over w,w, we obtain 



|Eq. {84a,)\'^ dw dw 



V^PlPl[vo X V2]2 



dtidiie'2''-*i~27„ti X 



X / dtodt'e'^(^'°~*'<^'>. 



to Uio 



(89) 



o 



(85d) We pass to the polar coordinates in the (ti,ti) plane: 

ti — tcoscj), ti — tsin(j), (90) 



Imr] ^ wvpiti+wvpiti + Ato, (86a) and parametrize the region O as 

A = (JWqi + 5ujq^ - SuJq^ - 6uJc,2 + 2Ae/i, (86b) 

O^{to: tCrmn{(i>) < to < tCmaxi^)} , (91) 

where all detunings at ti,ti have been absorbed into a 

shift of w,w. The energy mismatch due to electron-hole where Cmin{4') and C>naa;(</') are piecewise functions of 
asymmetry, Aeh = ao(pf„ — pLt); is analogous to that the form acostj) + PsiiKJ), corresponding to various con- 
defined in Eq. (74b). The damping factor is (we use the ditions (88). Performing the integration, we obtain 
fact that 70 = 70, 72 = 72): 

Re 77 = 70 (io + io) +71^1 +71*1 +72(^2 +^2) = ^ tt^ 



it/2 



fW2 X vi]^ , [vo X Vl]^ 

' 70 + 71 + T r 72 Ml + 



\[vo X V2]^ 
V2 X Vi], 

Vo X V2], 

Ixtl +Jytl- 



■ 70 + 71 + 



[V2 X Vo] 

[vo X vi]z 

[V2 X Vo]z 



^PlPlb^O X V2]2 A4 



72 t 



•. T\ /■ (M\ /■ (A\ 27a; COS </' + 27y sin (^ 

X T\ QmaxVP) - Cmm(0), T (P^a) 



(87) 



The two times ii,iican be taken to vary independently 
from to 00. The integration domain for to, which we 



T{x, y)= / (1 - cos xt)e~yh dt = ^^^^^i%^2b) 



J 
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< (/?2 < TT 



-7T < (^2 < 




TT TT + (fi2 27r 



TT + (/?2 TT 



FIG. 17: Region of integration over (pi,ipi, as determined by inequalities (95) for < ip2 < n and n < ip2 < 2ii. 

D. Angular integration and polarization dependence 

Let us collect all the factors in the expression for I2D* [we denote eo,2 = (cos (^0.2, sin 990,2)]: 

r |[eo X e,n]z? |[e2 x eout\z\^ x 

27a; COS + 27j, sin q 



dljK _ g /e2\ v^(\k\ ujl^t(ujin + t^outY f difio dip2 difii difii 



don 



2n 



\sm{ipo-ip2)\ 

7r/2 



. 2f0-fl .2'fil-f2 . 2^0 -'Pi . 2PI- V'2 

X sm sm sm sm 

2 2 2 2 

A system of two linear equations 



dcf) 

A 



--^^[C,nax{(f)~C,inin{4'),- ''^ "^ ' '« 



A 



Votm + Viti + V2to«i + Vifi = 



(,93) 



(94) 



has solutions with timtout,ti,ii > if and only if the four vectors vq, vi, V2, vi do not lie in one half-plane. Having 
found tin, tout, ti,ii, we can always split tin = to + io, tout =^2+^2 in such a way that to + ii + ^2 = ^0 + ^2-1- ^2 (since 
any side of a quadrangle is shorter than the sum of other three sides). This condition on the vectors Vo,Vi,V2,Vi 
translates into the following condition on the angles (since all conditions arc on ip2 — fa, Pi — fa. Pi — fa, we set 
(/?o = for brevity): 

< (^2 < TT : 

fi- [0;<P2] [P2;i^) {■k;p2+it) ((/32+7r;27r] 

Cpi: (7r;932+7i") (tt; ^1 + tt) (0; 27r) (931 - tt; (^2 + tt) 

— TT < (^2 < : 

fi : [0; ip2 + tt) {f2 + tt; tt] (tt; (^2 + 27r] [(^2 + 27r; 27r] 
(^1 : (<P2 + tt; (^1 + tt) (0; 27r) ((^i — tt; tt) {ip2+'n;Tr) 
<P2 = TT : 

(^1: (0;^) (7r;2^) 
(^1: (7r;2^) (0;7r) 

which define the domain of the angular integration. We also show it schematically on Fig. 17. 
Upon integration over (pi, Cpi Eq. (93) can be written in the form 



(95) 



47r 



dljK 

dOout 



dipodip2J{fa - P2) |[eo x ein]z\ |[e2 x eout]z\ ■ 



(96) 
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The function J(^p) is real and even, J[^^ = J{—^p), as can be seen from Eq. (78), so its Fourier series reads as 



27r 



J{^) = T, ^y^ -^ cos nip, Jn^ dip J (if) COS rvp. (97) 

ZTT ^ — ' TT J 



Then the angular integration gives 



^iK „_|_ |2| |2j»7o i/2r 2/ \ -2/ 



47r-- = 27r|ei„| \eout\ "^-r + "^[cos {p in - ^ out) - sin {p in ~ font)] 

aoout L * ° 

= 2^ (^|l - ^) \e,n?Kut\' + 2^ ^ [|e,„ • e:„,|2 + |e^^ . e„„,|2] . (98) 

Comparing this with Eq. (43), we obtain I^*^' = 27r(Jo/4 T Ji/^)- 

E. Dispersionless phonons 

For A = (electron-hole symmetric case) the last integral in Eq. (93) can be written as 

, 3 [Cmaxi4>) - CminW] 
9 4- 

16 [7a; cos (/) + 7j, sin 0] 



(99) 



Numerical evaluation of the angular integral (for simplicity we set all 7o = 7i = 7i = 72 = 7) gives the following 
expression for j7o, J2 to be substituted in Eq. (98): 

Jo\ _ fe^V^^ fXKYi^luMn+^out? 3^ / 0.0440, ,^^^ 



J2 J \c J c^ \2tt J 4 874 j^ -0.0017. 

This results in I-^j^ w ^Ik' ^^ ^^'^ polarization memory is almost completely lost. Adding the contributions from two 
mutually perpendicular detection polarizations [Eq. (48)], we obtain 

Uk = 0.0440 f^V 4 (—) ' "°"^"" ^ "°"^' J^ 4 - 3 cos 6,^. - cos^ 6,., ^^^^^,^ ^^^^^ 



2tt J 4 874 

For |ei„| ~ 1, Qdet = tt, uJin ~ ojout this expression corresponds^^ to Eq. (9) of Ref. 28. 

For |Aeh| 3> Jx,ly (strong electron-hole asymmetry) the last integral in Eq. (93) can be written as 

7r/2 



4A2 [7^. cos <j> + jy sin (f)] 



(102) 



Numerical evaluation of the angular integral gives 



'^o\_f(^^\ V^ f>^KY ^luti^in+^outf TT J2.60, 



J2 j \c J c^ \2n J 4 872 A^^ I 0.06. 

To obtain the expression for Uk one should replace 0.0440 —> 2.60, 7r2/(87'^) —> tt^ / {24j^ A^,^) in Eq. (101) 

e^X v'^f\K\ Ujl^t{uJin+i^outY ^^ 4 - 3 cos Odet - COS^ Odet I ,2 



(103) 



/4K = 2.60(-) -(^j °"""V """ 2472^2/ 4 ^"^"- |e.„|-. (104) 



F. Dispersive phonons manifold, the appropriate way to approximate the last 

In the case when A varies strongly compared to 7, 
but vanishes on a certain submanifold of the resonant 
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integral in Eq. (93) is 



7r/2 

27r(5(A) f a 



) [7a;COS(/) + 7j;Sin0]'' 



(105) 



This expression corresponds to complete suppression of 
interference between trajectories of different shape by 
the phase mismatch coming from the difference of the 
phonon frequencies. In this case the electron-hole dy- 
namics can be described by a kinetic equation, analyzed 
in Appendix C (up to an overall interference factor, see 
the discussion in the end of Appendix C) . 

First, let us focus on the apparent singularity at 
sin|(p2 — </5o| — > in the angular integral in Eq. (93). 
For generic ipi,(pi^ we see from Eq. (87) that ^^.y oc 
l/sin(v32 - </5o), and from Eqs. (88) - that Cmaxi4>) 
either stays finite (at ip2 ~ </'o), or also diverges as 
l/sin((p2 — 'Po) (at 9^2 ~ V'o + "■)• Thus, the power of 
sin((p2 — </5o) in the numerator of expression (105) is suf- 
ficient to suppress the singularity. 

The real danger in the angular integral comes from sin- 
gularities of the Jacobian corresponding to the resolution 
of (5(A), i. e., values of (pi, (f>\ such that 



^(Wpi 



Pi— Po 



Looking at the integrals 



'P2-Pl j _ rj ^l^Pl-po I ^P2-piJ _ n 



dpi 



(106) 
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FIG. 18: (color online). Plot of the left-hand side of 
Eq. (108a) as a function of ipi for tpn = Q, (p2 = 0.2 (a) and 
ip2 = T^ — 0.2 (b). The incident laser frequency ujin — 2 eV, 
ujo = 0.17 eV. 



dx dy S{x — 2/ ) = 



dx 

2xAx 



-2, 



(107a) Around each singular point we can expand 



27rrdr^==-,(107b) 



/ dxdy5{x^ + y^) = 

-1 

1 1 

jdxdy5{x'~y') = 2J^^=2\u^-, (107c) 



we notice that the logarithmic divergence appears when 
both ipi and (pi lie near one of these special points. Note 
that both ifi and (pi should lie near the same solution 
for the energy (5-function itself to be satisfied. 

At this point we have to assume a particular form of 
the phonon dispersion Wq. We take the conical dispersion 
cjq = 0^0 + Vp/i|q| valid for ujq/v ^ g <C l/a.^" Then the 
singular points are determined from the equation 

— sin((^i — Lpo) -\ sin((/?i — (^2) = 0, (108a) 

qi 92 



91 



92 



yPo + -Pi - 2poPi cos((^i - .^0), (108b) 
Y Pi +P2- 2pip2 cos(932 - fi)- (108c) 



which is the same for (pi and pi, since pi — pi- Let us 
denote the singular points by pi, labelled by the index s. 



^qi + ^^2 ~ ^s + ^ pI, (^1 = (/3i -(/?J,(109a) 

^qi + ^^^2 ~ ^s + -^ <pI, pi=<pi- pl.{109h) 

The second derivative in each singular point is given by 

-^ PiPo cos(935; - V3o) + -^ P1P2 cos(.^i -P2)- 
Zqi Zq2 



"^p'l 2 2 • 2/ s N "^ph 2 2-2/ 
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(110) 



To get an idea of the location of singular points, we plot 
the left-hand side of Eq. (108a) as a function of ipi for 
P2 close to pq and to (/Sq + "' (Fig. 18). We see that in 
the first case Eq. (108a) has two solutions, while in the 
second case - four. This means that there is a special 
value of (^2 between and tt (together with the sym- 
metric one), such that Eq. (108a) has three solutions 
and at the third solution the derivative w", defined in 
Eq. (110), vanishes. This situation will take place as long 
as LOo/uJin < (3 - y5)/4 = 0.191 . . ., i. e., uj^n > 0.9 eV, 
as can be established by setting ip^ = 0, pi — p2 = t^ [so 
that Eq. (108a) is satisfied], and equating uj'J — 0. How- 
ever, when these conditions are fulfilled, and when ip2 
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takes this special value, the third solution for tpi always 
lies in the smaller sector between (pQ and ip2, and does 
not belong to the integration region. As a result, there 
is always just one solution of Eq. (108a), which satisfies 

TT < (fii — ipQ < TT + ip2 — ifio ior < (f2 — fa < T^ and 

(^2 — <P0 — TT < (/5l — <po < TT for TT < </?2 — <P0 < StT [1. C, 

inequalities (95)], and thus contributes to the intensity 
integral. We denote this solution by (pi. 

Integration over the deviations (pi — pi 
'fi - fl gives 



dipi dipi 27r(5(A) — 
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Generally, the upper and lower integration limits here 
are of the order of ±1; more precise knowledge is not 
needed for for the calculation of the leading logarithmic 
term. It is important that the energy (5-function has a 
finite width ^ 7, which may cut off the divergency first, 
if it is greater than l^eh- The innermost integral over 0, 
assumed to be a non-singular function of (^1, (pi, can be 
taken at (pi = (^1 ~ p>\- To check the validity of this 
assumption, we have to study in more detail the behavior 
of the integral at p}2 — Po ^ 0, it. 

To make the formulas more compact, in the following 
we set Lpa = 0, as all angles can be counted from </?o- For 
simplicity we also perform the calculations in the limit 
ujq <^ LOin. Then, assuming — tt < p2 < tt, we simply 
obtain (pi — tt + P2l^- 

Let us start from the simpler case of (y92 close to ±7r, 
denoting (^2 = </52 + tt if — tt < (/32 < and <p>2 — '~P2 — t^ 
if < <y92 < TT. Then for <p\ = (p\ = ipo we have 
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2sin|(^2/2| 
H — min{— cosc/i + siiKp , 0} , (112a) 

Cmin{(f>) == -max{-cos(/) + sin0, 0}. (112b) 

For \ip2\ < 1 the condition Cmax{<P) > Cmm(0) severely 
restricts the integration domain in 0, so (p2 ~ ±7r does 
not introduce any extra singularities. 
For \(p2\ -^ 1 we have 
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Again setting all 70 = 71 = 71 = 72 = 7, we obtain sim- 
ply lx= ly = 27. The requirement C,min{<l>) < Cmax{(f>) 
translates into 



\hm < 



cost 



■ smc 



(114) 



If |(pi|,|<^i| ^ |</?2|, then only (pi « —ipi contribute to 
integral (111) (for tpi « ip>i the domain of integration is 
restricted by \(p>i\ < (/?2/2). Then constraint (114) al- 
lows only small deviations of (j) from 7r/4, thus we can 
approximate fi{(j)) ~ —^/2{pi/(p2){4'— 7r/4). In the op- 
posite limiting case, |(pi|, |(^i| <C |(p2|, range of 4> is al- 
most unrestricted, and both (fi « ±(^1 will contribute. 
The 0-integral in Eq. (105) in these two limiting cases is 
calculated to be 
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(115) 
Thus, we conclude that for \(p2\ <C 1 the upper cutoff 
in the logaritmic integral (111) is 1(^2 1, and not of the 
order of 1. Thus, upon integration over (p2 we obtain the 
second logarithmic divergence, which should be cut off at 
\(p2\ ^ 1 above and \(p2\ ^ max{7, Ae/il/fx"" below. 

We restrict ourselves to the calculation of the leading 
logarithmic asymptotics, so the second derivative (110) 
can be taken at (p2 ~ po = Q and is simply uj" = 
—uJinVph/{Av). The function J{p), defined in Eq. (96), 
can be taken to be 



J{V) 



^2 

8 


V 


e2\ 


c2 \2tt ) 


in 

1' 


V 
Vph 


X 


1 


-In 


ML^tnVph/v 

max{7, Ae/i} 







(116) 



The final integration leads to the following expression for 
Jo 7 1/2 to be substituted in Eq. (98): 
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so that IIj^ — 3/^, as for two-phonon scattering. Thus, 
the polarization dependence is the same as that described 
in Sec. VIA. Using Eq. (47), we obtain the final result: 
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To conclude this section, we note that the leading loga- 
rithmic term, calculated here, is not sensitive to the as- 
sumption of the conical phonon dispersion i^q. The same 
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FIG. 19: Electronic sell-energy Ironi the RPA-screened 
Coulomb interaction. 



result will be obtained for dispersion of any shape; the 
phonon group velocity should be taken at the wave vector 
q =^ ujin/v, corresponding to electron backscattering. 



VIII. RENORMALIZATION OF THE 
COUPLING CONSTANTS 



Coulomb renormalization 



As discussed in the beginning of Sec. VI B, the mea- 
sured ratio of the integrated intensities of 2K and 2r 
peaks at 2700 cm-^ and 3250 cm~i, hx/hr ~ 20 
(Ref. 6) is in noticeable disagreement with the calculated 
values of electron-phonon coupling constants. ^^ In this 
section we investigate how the electron-phonon coupling 
constants are renormalized by the Coulomb interaction 
between electrons. A brief account of this part has been 
reported in the short publication (Ref. 39). 

We are going to consider only the long-range part of the 
Coulomb interaction (i. e., smooth on the length scale of 
the lattice constant). Such interaction docs not mix the 
states in different valleys, so the interaction hamiltonian 
can be written as 



H,, 



-ld^v<fv'P^^P^''\ 



p{v) =il)\v)-ip{v). 



(119) 
In this equation we have not included explicitly the 
screening by the background dielectric constant of the 
substrate Eoo (the high-frequency value), which can be 
taken into account by incorporating it into e^. 

The electronic self-energy due to the Coulomb interac- 
tion (the Fock term), S'^'^(p, e) is shown in Fig. 19. Its 



leading logarithmic asymptotics is given by 
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S-(p,e)=^|^^y(q,^)G(p-q,e-.) 
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(120) 
arccosq ~, , q arccos q , 

, /(.g) = n — 1 ,(121) 
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(122) 



Here we have introduced the total number of the Dirac 
species, A/" = 4, which takes into account the valley and 
the spin degeneracy (the latter enters as factor of 2 mul- 
tiplying the electron polarization operator). The lower 
cutoff ^mm ~ max{wp, e}, the upper cutoff ^max ~ vja is 
of the order of the electronic bandwidth, and the dimen- 
sionless Coulomb coupling constant is defined as 
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8v 



(123) 



The derivation of Eq. (120) is given in Appendix D. The 
logarithmic divergence in the Fock self-energy T,^'^ is due 
to the long-distance nature of the Coulomb interaction, 
and thus is not picked up by local approximations such 
as LDA or GGA. 

The random phase approximation (RPA) for y(q, cj), 
shown in Fig. 19, corresponds to expansion of the coef- 
ficient in front of the logarithm to the leading order in 
the parameter 1/Af — 0.25, assumed to be small. This 
assumption is better justified than the expansion in g, 
which would be obtained if we used the bare coupling 
27re^/(7 instead of the RPA-dressed one V{q,uj). Indeed, 
for A/" = 4 we have g = (7r/2)(e^/w) w 3.4; taking into 
account the background dielectric screening reduces it to 

The presence of the large logarithm invalidates 
the simple first-order expansion in 1/Af, and makes 
it necessary to sum all leading logaritmic terms ~ 
{l/Af)^\n"{S,max/^min) of the pcrturbation theory. This 
summation is performed using the standard renormaliza- 
tion group (RG) procedure. Let us introduce the running 
cutoff ^maajC^^. Onc RG Step consists of reducing the cut- 
off, ^ ^ £ -h 6i, so that e*'^ > 1, while {1/Af)6e < 1. The 
inverse Green's function transforms as 



e-«p.S-I](p,e) 



e — (v + 5v)p ■ S 
1 + 6Z 



(124) 



where Z is chosen to preserve the coefficient at e upon 
rescaling of the electronic fields, ip ^ {1 + SZ/2)ip: 



= 1 
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'd^' 



1 + SZ 
The renormalization of the velocity is then given by 

Sv _ dY, as 

V de d{v\i ■ S) ' 
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FIG. 20: Diagram describing the logarithmic correction to the 
electron-phonon vertex due to the Coulomb interaction in the 
order 0{l/J\f). 
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FIG. 21: Plot of the function /(g), defined in Eq. (121) 



FIG. 22: Dependence of the dimensionless coupling constant 
\k on the running cutoff (,maxe~^ (to be identified with the 
electronic energy), represented by three upper curves starting 
from the bare value 0.04 at 10 eV, for three values of the bare 
Coulomb coupling g{Q) = 3.4, 1.5, 0.5 (solid, dashed, and dot- 
ted curves, corresponding to the substrate dielectric constant 
Eoo = 1, 2.3, 6.8, respectively) as determined by Eqs. (128). 
The constant Ar = 0.035 is unchanged. 



Next, we need to determine renormalization of the cou- 
pling constants. The electron charge is not renornializcd, 
as guaranteed by the gauge invariance, so the renormal- 
ization of the Coulomb coupling constant g is determined 
by renormalization of the velocity v. The correction to 
the electron-phonon coupling constants is determined by 
the diagram in Fig. 20, and is evaluated in Appendix D. 
The result is 
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(127b) 



Let us pass to dimensionless electron-phonon coupling 
constants Ar, A^, introduced in Eq. (24). Then the equa- 
tions for the RG flow are the following: 
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As fig) is positive and monotonous (see Fig. 21), 
g flows to weak coupling:"^^ if the initial value of g is 
large, 

2,9 

(129) 
while at small g we have 
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Integration of Eqs. (128) yields the following relation: 
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which, in principle, can be checked experimentally. Thus, 
Xk is enhanced, which is in qualitative agreement with 
the Raman data: according to the results of Sec. VIA, 
the ratio of the intensities of the two-phonon peaks is 
l2K/l2r^2{XAjXr)^ 

To study the behavior of the coupling constants quan- 
titatively, we solve Eqs. (128) numerically. The largest 
value of i is determined by the lower cutoff S,min ^ 
u!f^ ^ 0.2 eV. In Fig. 22, we show the flow of Xk 
for three values of the bare Coulomb coupling constant: 
5(0) — 3.4 (corresponding to no dielectric screening at 
all), g{0) = 1.5, and g(Q) — 0.5. The bare values of the 
the electron-phonon coupling constants Ar(0) ~ 0.035, 
Ak(0) = 0.040 were chosen (i) to satisfy the relation 
Ar(0)/Aif(0) — ujk/i^t, valid in the tight-binding ap- 
proximation, (ii) to reproduce the experimental value 
Ar ~ 0.035. Note that the RPA calculation without the 
RG collection of all leading logarithmic terms, would give 
all dependencies on Fig. 22 to be straight lines with slopes 
fixed at 10 eV. A comparable error would be produced 
by the GW approximation, which neglects vertex correc- 
tions, and thus picks up correctly only the first term of 
the logarithmic series. 

To estimate the EPC strength relevant for Raman scat- 
tering, we identify the running cutoff with the typical 
electronic energy, involved in the process, thus stopping 
the RG flow at electronic energies e ~ 1 eV (half of 
the incident laser frequency). In the unscreened case, 
(7(0) = 3.4, it gives Xa^/Xy ~ 3.2, in agreement with the 
observed ratio I2K/I2T ~ 20. 



36 



Finally, we wish to note that the cancellation of the 
self-energy and vertex corrections, leading to dXr/di = 
in Eq. (128b), is not occasional. Indeed, comparing 
Eqs. (23) and (28), we can see that coupling to the 
E2 phonon displacement u^;^ has the same form as the 
coupling to the vector potential A with the correspon- 
dence UE2X <^ Ay, —UE2V ^^ ^x (^P ^° ^^^ sign, differ- 
ent for the K,K' valleys). This means that a uniform 
phonon displacement u^;^ can be gauged out of the elec- 
tronic hamiltonian, which should hold for both intial and 
renormalized hamiltonians. Thus, gauge invariance re- 
quires that Fr is renormalized in the same way as the 
velocity v. Since Ar oc _Fp/u^, it must remain constant. 



B. Renormalization due to the electron-phonon 
coupling 

It turns out that the Coulomb interaction is not the 
only source of renormalizations. Since the electron- 
phonon self-energy ^^''(p, e), calculated in Sec. IV C, also 
has a logarithmic divergence, renormalizations due to 
electron-phonon interaction should be taken into account 
as well. However, in practice, the electron-phonon cou- 
pling is so weak (A^ <C 1, see the previous subsection), 
that its effect is negligible, so that the previous subsec- 
tion contains all the practical information. Still, for the 
sake of completeness, in this subsection we describe the 
theory of renormalizations due to electron-phonon cou- 
pling. 

Let us return to the electron-phonon self-energy 
SP''(p,e) = i;^(p,e) -hS^(p,e), calculated in Sec. IV C. 
The leading logarithmic asymptotics of T,^'^ is given by 
(see also Appendix D): 
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Here the phonon mode index ^ runs over the two modes 
belonging to the E2 representation and the two modes 
belonging to the E'-^ representation, the corresponding 
matrices (AE)^ being -A^Ej^, A^E^, A^-E^, A^^E^. The 
upper and lower cutoff is given by S,max ~ "^/a, £,min ~ 
max{e, w^}. The dimensionless constants Ar, A^, defined 
in Eq. (24), will be treated as small parameters. 

The latter statement deserves some discussion. In prin- 
ciple, one could proceed analogously to the Coulomb 
case: instead of doing the pcrturbative expansion in A^, 
one could dress the bare phonon propagators by the ap- 
propriate polarization operators n(q, w), corresponding 
to l/M expansion (the polarization operators for differ- 
ent matrix vertices are calculated in Appendix E) . Since 
n(q, a;) oc <7 at large q, the dressed phonon frequency 
would grow as ,Jq, and TP^ would no longer diverge 
logarithmically. However, the inelastic X-ray scattering 




FIG. 23: Diagrams describing the logarithmic correction to 
the electron-phonon vertex due to the electron-phonon inter- 
action in the order 0(1/7V, A^). Diagrams (b)-(e) vanish. 





FIG. 24: Logarithmic diagrams of the order O(A^) not re- 
duced to a renormahzation of the electron-phonon vertex. 



data for the phonon dispersion"^*^ show that the phonon 
dispersion is smaller than the phonon frequency itself. 
Thus, the renormalization of the phonon frequency re- 
mains small even at g ~ 1/a, so the perturbative expan- 
sion in A^ is more justified, and we neglect the phonon 
dispersion. 

The logarithmically divergent integrals in Eqs. (120) 
and (132) have different structure due to different form 
of the screened interaction I^(q, uj) and the phonon prop- 
agator D^{u)). In Eq. (120) the integral is dominated by 
the frequencies \lo\ ~ vq, while in Eq. (132) it is \uj\ ^ to^, 
since Dfj_{uj) ex 1/uj'^ at \uj\ 3> tj^. Thus, in the calcula- 
tion of the leading logarithmic asymptotics it is sufficient 
to approximate D^{uj) w — 27ri(5(a;) (see Appendix D). 
This substitution makes the phonon propagator (com- 
bined with electron-phonon vertices) formally analogous 
to the correlator of a static disorder potential (in other 
words, from the point of view of electrons with e 3> w^ 
the lattice is effectively frozen). Thus, renormalizations 
due to electron-phonon interaction at e 3> tj^ are equiv- 
alent to those due to static disorder^^"^^'^^. This equiva- 
lence holds only in the leading order in electron-phonon 
coupling, since in higher orders the phonon propagator 
is dressed by polarization loops, and the static disorder 
correlator is not. 

Logarithmic corrections to the electron-phonon vertex 
in the order 0{1/J\f, A^) are shown in Fig. 23. All dia- 
grams vanish, except the first one, which gives a nonzero 
correction to Xk' 
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The diagrams of Fig. 23, however, do not exhaust all 
logarithmic vertex corrections. In addition, one has to 
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consider two diagrams, shown in Fig. 24, as they are of 
the same order O(A^), and also logarithmically divergent. 
They may be viewed as a correction to the two-electron 
vertex r'^^)(p, e,p', e'; q, cj). The bare value of F^^^ is 



given just by the phonon single phonon propagator, com- 
bined with the electron-phonon vertices. Diagrams of 
Fig. 24 are evaluated in Appendix D, and give: 
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The following features of the expression for (JF^^^ are equations (128) as follows: 

worth noting: (i) the matrix structure of ST^'^' is identical 

to that of F^^); (ii) ST^^^ depends on electronic energies 

and momenta only through ^min, i- e., logarithmically; __E. — ■' ^^' 

(iii) the w dependence of JF^^^ is analogous to that of F'^^^ g M n'^Af 

but the pole is at the sum of two phonon frequencies. 



+ 2^ — ^ , (135a) 



n,k=0 



The latter fact has a simple physical meaning. Be- 
fore the reduction of the ultraviolet cutoff (,max the 
excitations of the lattice were phonons with momenta 
|q| < S.max/'v, as well as their combinations. In the the- 
ory with the reduced cutoff ^maxe~^, besides phonons 
with momenta |q| < {^max/v)e"^, one has to consider 
also pairs of phonons with large and almost opposite mo- 
menta q' and q — q', so that the total momentum of the 
pair is |q| < £,maxe^^- Each of the phonons constituting 
the pair has {^max/vje^^ < |q'|, |q - q'| < ^max/v and 
thus has been integrated out. The pair, however, hav- 
ing small total momentum, has to be included into the 
low-energy theory as a single excitation. Thus, D^+^' {uj) 
has the meaning of the propagator of this excitation, and 
comparing expressions (134a) and (134b), one can define 
electron-two-phonon vertex by analogy with the electron- 
phonon one. 



Iterations of the RG procedure will generate electron 
coupling to excitations with larger number of phonons, 
hence all these excitations with frequencies nuY + n'uj'^ 
have to be included in the low-energy theory separately. 
Obviously, electron coupling to excitations of the type 
nV + 2kEi will have the same matrix structure as cou- 
pling to the E2 phonons, while coupling to excitations 
nV + {2k + l)Ei' - the same as to the E[ phonons. De- 
noting the dimensionless coupling constant for the exci- 
tation nujr + n'uj'^ by A„^„', we can generalize the RG 
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These equations can be simply related to Eqs. (12) of 
Rcf. 39, where Ar and A^i denoted the total oscillator 
strength of all E2-\ike and all E[-\ike excitations, respec- 
tively: 
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Namely, Eqs. (12b), (12c) of Ref. 39 are obtained by sum- 
ming Eqs. (135b), (135c) over n and k. 

We will not solve Eqs. (135), since in graphenc the ef- 
fect of the electron-phonon interaction turns out to be 
neghgibly small due to the smallness of Ar,A_R-, as we 
mentioned in the beginning of this subsection. In partic- 
ular, the modification of the plot in Fig. 22 would not be 
noticeable by naked eye. 
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IX. SUMMARY 

In this paper we have calculated the frequency- 
integrated intensities of two- and four-phonon Raman 
peaks in disorder-free graphene. We started by writ- 
ing down the low-energy hamiltonian of the interaction 
of electrons with the crystal vibrations and the electro- 
magnetic field from pure symmetry considerations; as a 
result, we obtained a description of the system in terms 
of just a few independent coupling constants, considered 
to be parameters of the theory. Another parameter of 
the theory, introduced phcnomcnologically, is the elec- 
tron scattering rate 27. 

First, we analyzed the onc-phonon peak at 1580 cm^^ 
and have shown that the scattering is completely off- 
resonant: the intermediate electron and hole states in the 
whole first Brillouin zone contribute to the Raman scat- 
tering amplitude. As a result, the intensity of the peak is 
expected to be insensitive to most external parameters: 
polarization, electron concentration, degree of disorder, 
etc. However, according to our results, it is proportional 
to the fourth power of the excitation frequency. 

Then we calculated the intensities I2K and /2r of the 
two-phonon peaks at 2700 cm"-'^ and 3250 cm~^, respec- 
tively. We have shown that two-phonon scattering is 
fully resonant, so that the intermediate electron and hole 
states correspond to real particles, propagating along the 
quasiclassical trajectories and subject to scattering pro- 
cesses. As a result, the intensities are determined by the 
electron scattering rate. Besides, the quasiclassical char- 
acter of the process imposes a severe restriction on the 
electron and hole trajectories which can contribute to the 
two-phonon Raman scattering: upon the phonon emis- 
sion the electron and the hole must be scattered back- 
wards. This restriction results in a significant polariza- 
tion memory: it is almost three times more probable for 
the scattered photon to have the same polarization as 
the incidcnd photon than to have the orthogonal polar- 
ization. 

We have also calculated the intensity I4K of the most 
intense four-phonon peak at 5400cm~^. The four-phonon 
Raman scattering is also fully resonant. As a conse- 
quence, we have shown that measurement of the ratio 
I4K/I2K enables one to extract information about the 
relative contributions of different processes to the elec- 
tron scttering rate. 

Having compared the experimental two-phonon peak 
intensities, we extracted the ratio of the corresponding 
electron-phonon coupling constants. This ratio turned 
out to be significantly different from that obtained earlier 
from the density-functional theory calculations. We have 
shown that the reason for this discrepancy is the renor- 
malization of the coupling constants due to the Coulomb 
interaction between electrons, missed by DFT calcula- 
tions based on local or semi-local approximations for the 
exchange-correlation functional. In particular, we found 
that the constant, responsible for the peak at 3250 cm~^, 
is enhanced, and this enhancement is in quantitative 



agreement with the experimental Raman data, provided 
that the screening of the Coulomb interaction by the sub- 
strate is weak. 
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APPENDIX A: MATRIX ALGEBRA 

The natural basis in the 16-dimensional space of all 
4x4 hermitian matrices is represented by the r/^^ x 
Tf^, I, J ^ x,y,z 



0, where t,^^' ,t^^, 



J 1 "1 J ' -^^ ^7 '^7 ^: " ^^^^^ • ^ 7 • J T '^7 J '^lU} ^7 are 
the Pauli matrices acting in the corresponding subspaces, 
respectively, and t^^ , t^^ are the 2x2 unit matrices. 
The square of each matrix is equal to 1 - the 4x4 unit 
matrix. 

These 16 matrices can be split into two sets: those di- 
agonal and those off-diagonal in the KK' subspace (i. e., 
containing Ti^^ jT^^ and t^^ ^'''y^ ■> respectively). 
These two sets are invariant with respect to Cgi, because 
any transformation from Cgu either (i) leaves K and K' 
in place, thus belonging to Cs^, or (ii) swaps between 
K and K' (belonging to C2Czv), but never mixes them. 
Thus, these two sets form two 8-dimensional representa- 
tions. Both are reduced as Ai + A2 + Bi + B2 + Ei + E2, 
since each set contains a matrix transforming according 
to Ai: for the valley-diagonal set it is 1, while for the off- 
diagonal one it is the matrix Uc2 of the C2 rotation which 
commutes with all Cq^ , and is hermitian since C2 = C2 ■ 

Let us focus on the valley-diagonal set. The two ma- 
trices which transform according to Ei (vector) rep- 
resentation will be denoted by Sa;,Sj,. To establish 
their multiplication rules, we form the direct product 
El X El = Ai + A2 + E2. The corresponding linear 
combinations are 



SojSic + ^v^v ^ Ai, 




(AI) 


I^x'^y — SyEa; ^ A2, 




(A2) 


{'Sx'^x — ^y'^y, Sa;Sy + 'EyT^x} ^ 


-E2. 


(A3) 



Since Yc, 



1, the proportionality coefficient in the 



last line must be zero, so Ti^^y = — SyS^, = iS^. The 
matrix S^ defined in this way (i) is hermitian, (ii) trans- 
forms according to A2, and (iii) S^ = 1. Thus, it must 
coincide with the corresponding matrix from the basis 
(the sign may need to be changed). In other words, the 
set {Sa;, Sy, S^} satisfies the usual Pauli matrix algebra. 
As the subgroup C2V = {E, cr„ , cr(, , C2 } is a direct prod- 
uct, the matrices of i?i,i?2,^2 representations of the 
valley-diagonal set must commute, and the product of 
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any two will give the third one. Let us denote by A^ the 
matrix of the Bi representation, then that of B2 is A^E^. 

The two matrices transforming according to the E2 
(pseudovector or tensor) representation are denoted by 
{T^,Ty} = T. Just like for S, using E2 x E2 = 
Ai + A2 + E2, we establish the Pauli matrix algebra 
for the set {Tx,Ty,T,z}- It is important that matrices 
of cr'a^'^bJ '^'c reflections are expressed in terms of T^^Ty. 
Indeed, (i) a'^ ~ (cr^)^^, so its matrix is hermitian and 
its square is equal to 1; (ii) this matrix is diagonal in 
the KK' indices; (iii) C^^a'JJ^j, = cr^, so the matrix 
of a'^ cannot belong to a one-dimensional representation; 
(iv) C2<7'aC2 = <y'ai SO it must bc the E2 representation. 
Since we have a rotational arbitrariness in the choice of 
Tx, Ty, we simply fix Ty to be the matrix of the cr^ reflec- 
tion. 

The products of S and T can be analyzed by using 
Bi X El = B2 X El — E2- The rotational arbitrari- 
ness in the choice of Y^r^^J^y can be removed by flxing 
their behavior under the a'^^ reflection: TyT^^Ty = S^,, 
TyYiyTy = ~Tj y . Thc algcbralc relations established ear- 
lier leave us with 



T, --CA.Sj,, T^^CA.S,, C =1 



(A4) 



The sign of C, is not important as one can always redeflne 
A2 -^ — A2. The difference between real and imaginary C 
is essential, as it determines the symmetry under the time 
reversal, and will be discussed below. 

The off-diagonal set can be obtained from the diag- 
onal one by simply multiplying it by C/c'27 the matrix 
of C2 , which transforms according to the identical repre- 
sentation of Cgt;, but swaps K and K' . According to thc 
representation algebra, Uc2 commutes with T,z,Tx,Ty, 
and anticommutes with A^, A^Sz, S^;, Sj,. Generally, if 
two hermitian operators commute, their product is her- 
mitian, while if they anticommute, their product is an- 
tihermitian, so thc proportionality coefficient in the cor- 
responding algebraic relation must be imaginary. Hence, 
the matrices of the off-diagonal set can be written as Uc2 7 
iUc2-^z, Uc2^z, iUc2^z^z, iUc2^x, iUc2^v, -Uc2Tx, 
—Uc2Ty Denoting [/pjE^ = A^:, we recover the bottom 
row of Table III, provided that C = 1- 

To establish the form of the matrix of the C3 rotation, 
we note that (i) it must be diagonal in the valley sub- 
space, (ii) commute with Uc2 ■, (ii) its third power should 
be equal to 1. This fixes Uc^ = g±(2^i/3)s,^ 

Now let us establish the symmetry properties of the 
matrices with respect to thc time reversal. The explicit 
form of thc corresponding matrix Ut need not be spec- 
ified. We only notice that time reversal must commute 
with any spatial transformation. Applying this condition 
to the matrices Uc2^ T, e"^^^^'/^, we obtain Uc2 ^^ Uc2, 
T I— > T, S^ 1-^ — Sz. Next, the time reversal swaps 
K and K', so A^ 1-^ — Az. Applying the algebraic rela- 
tions established above, we obtain 



A, 



r 



'A, 



Sz' 



Sh 



-C^s. 



(A5) 



Next, since the matrices of (Ja,(Jb,<^c reflections are ex- 
pressed in terms of C/C2T, we must have Aj,S h^ Aj,S, 
which fixes ^ = ±1. 

To help those readers who prefer to work with a par- 
ticular representation, rather than basis-independent al- 
gebraic relations, we give specific expressions for the ma- 
trices defined above. In Ref. 54 the representation is 
introduced by defining the column state vector as 



i> 



IpAK 
IpBK 
IpBK' 

IpAK' 



(A6) 



and thc Si , S' matrices have the form: 



s. 


-r^^ 


' X I 


A, 


- r^^'r^"" 


s. 


-rf^ 


^AB 

ly , 


A, 


^ ^KK' AB 

'y 'z 


Sz 


^r!f^ 


rf^, 


Az 


KK' AB 



(A7) 

where r/^^ 7 '''j 7 ^7 J = 2;,y, z, are the Pauli matrices 
acting in the corresponding subspaces, respectively, and 



^KK' 



Tn^ are the 2x2 unit matrices. 



In Ref. 53 the column state vector is defined as 



^ = 



IpBK 

IpBK' 

-IpAK' 



(A8) 



and it is assumed that Wk'.a(b)(>^) — ^k AiB)^^)' ^^ ^^^^ 
representation the electronic matrices acquire an espe- 



cially simple form: Ej 

the time reversal matrix Ur 



^KK'^AB 

'0 ~i 7 



A,; 



^KK'^AB 
^0 



and 



KK' AB 

y y 



APPENDIX B: EFFECTIVE HAMILTONIAN IN 
AN EXTERNAL FIELD 

Thc theory of electrons in a crystal lattice subject to 
a magnetic field was developed long ago.^^ Here we con- 
sider the specific case of the two-dimensional graphene 
crystal. 

Let us start from the simpler case of the scalar po- 
tential (^(r, z), assumed to be smooth on the scale of the 
lattice constant. The definition (9) of the effective hamil- 
tonian can be written as 



■Pli^)Hlfi-pfi{v)d^v= / ■pi{v)eLp{Y,z)iPp{Y) X 



A^,-^^)rui{v,z)Up{v,z) 



d vdz. 



(BI) 



Here a, (3 = 1,2,3,4 label the four states with zero en- 
ergy, Kq^^ are K or K', correspondingly, and summation 
over repeating indices is assumed hereafter, unless stated 
explicitly otherwise. 
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The integration over z in Eq. (Bl) is straightforward 
due to confinement provided by the Bloch functions. As 
for the integration over x, y, we note that the combi- 
nation of functions in the square brackets can be repre- 
sented as 8apXa^'^)~^Xa-(i(j^ -z), whcrc the second function 
(i) has zero spatial average, and (ii) is periodic in x, y 
with the period corresponding to the tripled unit cell. 
The first term is written using the normalization (7) of 
the Bloch functions (no summation over a is assumed). 
The rest of the integrand, ^*JyY) (/3(r, z) V'/3(r), is a smooth 
function of x, y, i. e. its spatial harmonics have wave vec- 
tors much smaller than the inverse lattice constant. 

To study the integral of an arbitrary smooth func- 
tion /(r) with the periodic function Xa/3(r, z), we expand 
the latter in the Fourier sum over the reciprocal lattice 
vectors G: 



X.Mr,^)=E E Cip{<2,,z)e<'^^^'^>. (B2) 
Then the integral can be rewritten in the Fourier space: 



/(r)x./3(r,z)d2r = )J Y. ,f(G + jK) C^^(G, z) 



(B3) 



E E 

G i=0,±l 



The sum does not contain the term with G = 0, j = 0, 
which has been excluded from Xafi{^^ z) by construction. 
The Fourier transform of /(r). 



/(P) 



-ipr 



f{r)d\, 



(B4) 



is rapidly decaying away from p = 0. If all spatial deriva- 
tives of /(r) are continuous, this decay is exponential. 
This leads to the following expression for the effective 
hamiltonian: 



ttV 



eip{r, z = 0)Sai3 - 4f^(r, z = 0)(5a/3,(B5a) 
tPrdz 



ez\Ua{Y,z)Y 



LxLy 



(B5b) 



where E^ = —dz(p. The fact that dz does not depend on a 
follows from the transformation properties of the Bloch 
functions under the cr t,,(T(, reflections. The matrix ele- 
ment between states with Kq, =^ K^ is small as e~^^^P°'\ 
and thus cannot be included in the regular expansion of 
the effective hamiltonian in the parameter pa <C 1 . 



J 



To describe the effect of an external magnetic field, we introduce the vector potential A(r, z) in the gauge A^ — 0. 
The microscopic hamiltonian is given by 



Ti"^ = -— [A,(r, z) Vi + V, A,{r, z)] 
Zc 

In the first order of the k • p perturbation theory we obtain 



2inc' 



■A\v,z), 



-iV 



(B6) 



-^ / ro.{v)Ul{v,z)e-'^-^[M{v,z)v,+v,M{v,z)]e'^^-Ufi{v,z)i^p{v)d^vdz 



= — {a\Vi 
c 



[a zvi 



d^rC(r)A,(r,0)^^(r)--r5„^ / rf^r CW [A,(r,0) «, + «,A,(r, 0) ] V/sW 
2 , 9A,(r,0) 



(B7) 



where we introduced the following notation for the matrix elements of an arbitrary operator O between the Bloch 
functions: 



{a\0\l3) 



,-iKo 



-U:{r,z)Oe''^^''Ufj{r,z) 



d^rdz 

LxJ-ly 



(B8) 



We again encounter integrals of smooth functions with periodic ones, such as U'^e^^^"'^ Vie^^'^'^ Up , so the matrix 
elements are different from zero only if Kq, = K^. Recalling that in the first order of the k • p perturbation theory 
the Dirac hamiltonian is obtained as 



[i/i(p)]a/3 = (a|v|/3) • P = vY^ap ■ P, 



(B9) 

-iV — (e/c)A in 



the first term on the right-hand side of Eq. (B7) can be identified with the gauge elongation p 
the Dirac hamiltonian, Eq. (28). 

In the second order of the k • p perturbation theory one has to include the contribution of remote bans b ^ a Va = 
1, 2, 3, 4 with energies Eh, taken at the points K, K'. The correction to the Dirac hamiltonian is then given by 



P X — ^' 

[-f^2(p)]a/3 = — 5ap+PiPj2_, 

b 



{a\v,\b){b\v^\P) + {i 
2Eb 



(BIO) 
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The gauge elongation p -^ — iV — (e/c)A of the first term on the right-hand side corresponds to the second term on 
the right-hand side of Eq. (B7) and the A^ term of Ti. . The effective mass tensor originating from remote bands in 
Eq. (BIO) is symmetrized with respect to i,j in order to ensure the hermiticity of the effective hamiltonian upon the 
gauge elongation. Taking Ti. in the first order (the contribution of the div A vanishes, as it is a smooth function), 
and PiVi in the first order, we obtain the following contribution of the remote bands to the effective hamiltonian: 

b 

The jj-symmetric part of this expression corresponds to the gauge elongation of the last term in Eq. (BIO), while the 
antisymmetric part can be rewritten as 

em ^^ ia^vWMm^llljl [^,(r, 0) v, + v, A,(r, 0)] = 



2c ^ Eb 



{a\v,\b){b\v,\P)-{z^j) 



17 2^ ' 



4c ^ Eb 

b 



dMr,0) dA,{r,0) 
dxj dxi 



(B12) 



The expression in the square brackets is an antisymmetric tensor whose xy and yx components are equal to Bz and 
—Bz, respectively. Thus, this term corresponds to —fizBz'^z term in Eq. (29). 

Since Bx = —dAy/dz, By = dA^/dz, the last term in Eq. (B7) corresponds to the ^xy{Bx^y — ByTix) term in 
Eq. (29). Using the facts that (i) {a\z\(i) = {dz/e)Safi, (ii) zvi — ViZ, (ui) (a|2|6) — {i/2Eb){a\vz\b) (the latter follows 
from the commutation relation [H, z] = —{i/2)vz, valid in the absence of the field), we can write 

{a\zvi\(i) = — {a\vi\(3) + 2 Z^ ^ = — (a|w»|p) - ^/^ e " ^ ' 

b b 

The contribution of the {dz/ e) {a\vi\[3) term to the effective hamiltonian has the form —dz{v'ExBy — vTiyBx) and may 
be viewed as the Lorentz force contribution to the —dz£z term. 

APPENDIX C: KINETIC EQUATION 

Kinetic equation is the most natural way to describe the dynamics when real quasiparticles (electrons and holes) 
move along quasiclassical trajectories, and interference between different trajectories is suppressed. Since electrons 
and holes are created in pairs, the kinetic equation should be written for the joint distribution function fp^.-ph ('"ei ^h] t) 
- the joint probability for the electron to be in the elementary volume of phase space d^re (i^Pe/(27r)^ around Tg, Pe 
and for the hole in the elementary volume d'^rhd'^ph/i^n)'^ around rh,Ph- We write the kinetic equation as 

'-'/ I C'J^Cc;p^ OJ O-ttCc;— p^ OJ Q.e.out r I Cfe,iri f , Qt-h.out r , aih,in r I r^^ n\ 

m + ^^a^ + ap, a^ ~ ^^^-^- ^ + ^^p-p" ^ + ^^^-^- ^ + ^^^-^- ^' ^^^^^ 

Stpe.Ph / = -^^ j J2^ J^ — sm 2 '^(I^c^Pe - ReCp.-q - ^q),/p.,p;,, (Clb) 

SVe,P. ^^^^ J J2^ J^ ~ 2 '^(^e^P=+q - ^^?Pe - ^q) /Pe+q,ph> (Clc) 



Kzi f = -^-l -^ iS;^ "°' ^"^ t"^" '^(^^^"-- - ^^^'-^^-^ "'^^ ^— ^^^^) 



P''-P'^-'~J72^aK) 4 — 2 '^(^c^-Ph-q-i^e^_p^ -tjq),/p^^p^+q. (Cle) 

The left-hand side of the kinetic equation (Liouville operator acting on the distribution function) represents the 
free propagation of the electron and the hole with the corresponding (group) velocities Vg = 9Re^p^/9pe and 
Vft, = dKe^-p^/dph- The right-hand side (collision integral) describes emission of phonons; assuming to be in 
the linear regime, we have neglected the Fermi statistics of electrons and holes. In what follows, we will write the 
out-scattering part of the collision integral as — 2(7p^ + J-p^)fp^,ph, and include it in the Liouville operator. In 
the scattering rates 27p^ , 27_p^ we also include those for emission of electron-hole pairs, or other excitations with 
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broad spectrum; the contribution of such processes to the in-scattering part of the coUision integral is neglected (see 
discussion in Sec. I A). 

The Raman signal is treated as a weak probe not affecting the electron and hole population, and thus is not included 
into the kinetic equation. It is determined by the radiative recombination rate, and can be calculated from the Fermi 
Golden Rule. The probability of emission of a photon with a given polarization Bout per unit solid angle, per unit 
frequency interval, and per unit time is expressed in terms of the joint distribution function as 

Audi fev\^2nc^ uj^^f. /" 2 (Pp 



, , ,, = 27r r^y^I^^ fd'r-^S{u^out - Rc^p - Rc^p) |[eo.* x ep],|Vp,-p(r,r). (C2) 

doout dujout dt \ c J ujout 2t:'^& J {2TTy 

The easiest way to arrive at this expression is to calculate the electronic radiative self-energy: 

Note that the imaginary shift of the poles in the electron Green's function is different from the prescription (31). 
Indeed, the latter corresponds to the full valence band and empty conduction band, while radiative recombination 
requires the valence band to be empty, hence the shift of both poles to the lower half-plane of e in Eq. (C3). 

Kinetic equation (CI) contains no term corresponding to the generation of electron- hole pairs by incoming photons. 
Instead, we prefer to choose the initial condition at t = 0, which would correspond to the electron-hole population 
upon arrival of a short pulse of the electromagnetic field of the total energy Wi„, i. e., containing a single photon: 

/p^.p,(re,r„;t = 0) = — 2|[e,„ x epj,]^ '^^'; ~ ""^^ ^ ^(Regp^ + Reg_p, - c.,„) (2^)^<5(pe + p^). (C4) 

The last (5-function takes care of momentum conservation during photon absorption: since the photon momentum 
is very small, the electron and the hole must have opposite momenta. The energy (5-function ensures that the total 
energy of the electron-hole pair is equal to the energy of the absorbed photon; the coefficient in front of it is just the 
inverse density of states of electron-hole pairs with zero total momentum, necessary to preserve the normalization of 
/pe,ph (''e, r/i). The factor 5{Ye — Vh)l{LxLy) reflects the fact that the electron and the hole are born at the same spatial 
point, which can be located anywhere in the sample. The angular factor 2| [ei„ x ep]^ p shows that the transition dipole 
is perpendicular to the electron momentum, and the factor of 2 fixes the average to unity. As a result, the integral of 
/p^.p^(re, r/i) over the whole phase space equals 7re^/c, which is nothing else but the total probability for the incident 
photon to be absorbed; all subsequent factors, discussed above, represent the partitioning of this probability over 
different states of the electron- hole pairs. The total probability can be found directly from the Fermi Golden Rule, 
and is given by 



v^E 



2 

27rc2 



VuJi, 



27r(5(Re^p-hReCp- w. 



in f ; 



where the factor of 4 keeps track of the valley and spin degeneracy, and Lz/ c is the attempt period. 

Solution of the kinetic equation (CI), corresponding to emission of n phonons, is obtained by n iterations of the 
collision integral. First of all, we find the inverse of the Liouville operator, acting on a source Jp^,p^(re, v^; t): 



Jp.,P,(re,r„;i) - / dt' J^^,^Jv,-v,(t-t'),Vh-^h{t-t'yA')e-^^'^'^''^'^^'~'"^ ■ 



| + 2(7e+7.)+Ve^+v,A 

— 00 

(C5) 
Let us follow evolution of a single electron-hole pair, created at i = at the point r = (evolution of the initial 
condition (C4) can be obtained by a simple convolution). The zero-approximation distribution function is given by: 

4°)p^(re, r„; t) - {2ti)H{p, - po) (27r)2,5(p„ + po) 6{v, - v,t) 5{vh - v^t) e-27po*-27po*. (C6) 

After one iteration of the in-scattering part of the collision integral we obtain the first-approximation correction - the 
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contribution from electrons and holes which have emitted one phonon: 



/«pjre,r;,;t) - f(^}^^^{r,,r^;t) + f^1]^^{r,,r^;t) , (C7a) 

4:)pjre, r„; t) - (27r)'^<5(p;, + po) <5(r;, + vot) 6-^7™* x 

sm^ PO p- ^(Rc^p^ ^ Rcgp_^ - a>p„_pj X 



^^Po-Pe 4 2 



X y dto <5(re - Ve(t - to) - Voto) e-27p.(*-to)-27poto^ (C7b) 



4^)pjre,r;,;t) = (27r)'^<5(pe - po) <5(re - voi) e-^^™* x 

^^ 27rf|- \/27^ 2 (^p„ - (^-p^ ^^^^^ - - 



^^Po-Pe 4 



X j dio 6{rh - vh{t - to) + voio) e-2^-p'>(*-*«)-2^po*"o. (C7c) 



After the second iteration we have^^ 

/(2)p^ (re, r;,; i) = f^l,^ (r„ r;,; t) + f^^^ (r^, r;,; t) + 4^1;)^ (r^, r;,; t) , (C8a) 

/^-),Jre,r;,;t) - (27r)'^<5(p;, + po) (5(r„ + voi) g-^^-* x 



/■ d^pi 27r^|. \/27a2 . 2 (^ - (^ e/o <- ^ <- n 

" i (2^0^ M^p„_p, 4 "" ^^ ^(^^^- - ^^^- - "---) " 

X JdhJ dto S{r, - Ve(t - h) - vi(<i - io) - voio)e-'^-'*"*^^"'^-'*^"*'''"'^™*", (C8b) 


t 
X /"dio'5(re-Ve(i-io)-voio)e-27Pe(*-*o)-27p„*o X 





27r^|. \/27^ 2 '/'po - '/'-Ph cm^ J p^ J , , N^ 

- d(Ke^p„ - Ke^_p^ - tj_p„_pj X 



^'^Po-Pe 4 2 



X y dio 5{Yh - ^h{t - to) + voto) e-2^-pJ*-*o)-27poto^ (C8c) 


/^^;;)^(re,r;,;t) = (27r)'^<5(pe - po) <5(re - vot) g-^Tpo* x 

/• d^pi 27rf|- V27^ . 2 y^po - ypi f,o F T5 e ^ 

^ y (2^ M^p„_p, 4 ^^" ^ ^(ReCpo -Roep. ---p„+pj x 

X ^ll'^_^^ ^ sin^ ^^^^^ ^(RoCp. - ReC-P. - -..-pj X 

X fdii /dfo<5(r?.- v,,(t-fi) + vi(fi -fo) + vofo)e-'^-'''^(*-*^'-'^pi(*^-*"^-'^''»*°.(C8d) 


After four iterations we have f^^^ = f^"""^ + /(^^^eft) _,_ j(eehh) _^ ^(e/ih/i) _^ j(/ih/ih)^ r^j^^ ^^^^ contributing to the 
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four-phonon Raman signal is given by 



f^:TJire,r,;t) = / |pL _M^^ sin^ ^^^^ ^iRc^p^ - Rc^p, - c.p„_pj x 



(2^)2 M^p„_p, 4 



27ri^?- v27a? o Wn, — (^n ^.„ ^ „ ^ 
X ^^J_^^ ^ sm^^^^^^(Roep. -Rc?p. -u.p,_pj X 





d^p^ 27ri^|. V27a2 . 2 "y^po - '/'pi rm ? r. ? \ 

" ' (2^ M^p„_p, 4 "" ^^'^(R^epo-ReC,,-c._p„+pJx 



2^F|- V27^ . 2 '/'pi - y-P. 



A^^-Pi-P. 4 



(5(RcCpi - RcC-p;. - w-pi-pj X 



X / dh / dto S{rh - Wh{t - h) + vi(fi - to) + voto)e-27-P,J*-*i)-27pi(ti-to)-27p„to, (C9) 


In terms of this correction, the Raman scattering probabihty can be expressed as 



47r- 



dl 



AK 



doo 



2tt 



2 27rc2 ujI„, 27re2 87rw2 



C / i^out 2tT^C^ C CiJ. 



X / d\ £^ S^\[e.n X epj,|2|[e™, x epJ^p^lRe^po + Re^"p„ - ^„0 4::'!p',(r,r). (CIO) 



(27r)2 (27r)2 
Let us change time integration variables according to 



dt 



t t-i 
dti I dtQT(t-ti,ti -to, to) 




t ti 
dti / dto Tit ~ti,ti-io,io) 




= / dtodti dt2dtodtidt2S{to + ti + t2 - to - ti - t2) T{t2,ti,to) T{t2,ti,to), 


and evaluate the time integral, following Sec. VII C. We obtain 

cx) 7r/2 



(Cll) 



— — i — — I dti dh e~27.*i-27,u r ^^^ ^ — L^ 

2|[voXV2j2|7 J |[voXV2jz|J 8[72,cos(/) + 7ySin0l' 

o " 

Rearranging the energy (5-functions and we arrive at 

4 



(C12) 
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dl4K 

dOrm.t 



X ei„l^|^|fv2 X eout]zf X 



^^^'(DX^) '''J^'Pod'pid'p,d^p2 |[voxe,„],|^|[ 

X sin^ ^i^^ sin^ ^1^^ sin^ ^" ' ^^ sin^ ^1^^ x 
2 _ 2 2 2 

X 2(5(Re^po +Refpo - uj„,)S (Re S_p^ -Rc^p^ + c^p^^pj (5(Rcfpj -Rcfp^ + w^p^+pj x 

X 2(5(ReCp2 +ReCp2 - ujin + (^po-pi +^pi-P2 +'^-po+pi +'^-pi+p2) x 

X '5(Re^p, -ReCp, - Re^po +ReCpo +'^po-pi +'^pi-p2 - ^^-po+pi -'^-pi+p2) x 

ir/2 

1 /" JJ, Cma£c(0) - Cmm(</') 
X 77 TT / "?' T • 

vo X V2 ^ 7 8 [j^ COS (/) + 7y sin 0] 



(C13) 



The first four (5-functions constrain the electronic momenta to lie on the resonant manifold (Fig. 16). The argument of 
the fifth (5- function is nothing but A, defined in Eq. (86b). For ujq = ujo this (5- function equals either 00 in electron- hole 
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q, uj 
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p,€ p + q 

e + LJ 



q,w 



p + q 

'7 


/ \ p + q + q' 


A 

p,e 


q,'-^ p + q 



FIG. 25: Coulomb correction to a generic vertex — ir(p, e; q, a;), shown by the triangle. 



symmetric case, or in the electron-hole asymmetric case. Thus, Eq. (C13) makes sense only for dispersive phonons. 
In this case it gives the probability 4 times smaller than that given by Eqs. (93), (105). One factor of 2 is due to the 
constructive interference of the amplitudes for the two spin projections: the spin degeneracy multiplies the matrix 
element, obtained by tracing in the electron loop. The other factor of two is due to the constructive interference of 
two time-reversed processes, described in the end of Sec. VII A, which may also be viewed as the interference of the 
processes in the two electronic valleys. 



APPENDIX D: LOGARITHMIC TERMS IN SELF-ENERGY AND VERTEX CORRECTIONS 



The RPA-dressed Coulomb propagator (122) determines the electron self-energy: 



„»». ^ f dcu (Pa , ^ , , /" dvj (Pa , , „, 

^'^P,e) ^ I J —-^Vi<i,u;)Gip - ci,e- cj) ^ - J —^ y(q, ztn) G(p - q, 6 - Jti7) 



16g [ dvo (P{vq) 1 \/{vqY + ^^ 



A/" 7 27r (27r)2 vq gvq + ^(wg)^ +■ 



ttW 



d^ / dzu 



ve 



9^ + ^e + tJ72 



27r (27r)2 
itu + TjqS e + f pS 2(ien7 — wpq)(in7 + i;qS) 



tu^ + {vq)'^ tu^ + [vq)"^ 



{w^ + («g)2)2 






wpS 



In- 






dip 



(-e + wpS) 



2(1 + gsimp) 



(2e - upS) 



g sm ip 



„, garccosg , ^, / 

i-e + vpi:) \ ^ + (2e - «pS) 1 - 

V 1 - r V 



2(1 + gsaiif) 
TT arccos g \ 



2g g^Jl^^ 



(Dl) 



In the first line we have performed the Wick rotation uj — izu (corresponding to the Matsubara representation for 
zero temperature). In the second line we have expanded G{p — q, e — iw) to the first order in e and vp, since the 
integral is dominated hy w ^ vq ^ e, vp. In the third line we have integrated over the directions of q. In the fourth 
line we have replaced the integration region with the momentum cutoff ^mm < vq < S.max, — oo < n7 < cxd by the 
region ^mm < ^ [vq)'^ + ^ < £,max, which does not change the leading logarithmic asymptotics. The lower cutoff 
Cmm ~ max{wp, e}, the upper cutoff S.max ^ v j a is of the order of the electronic bandwidth. 



Let us now calculate the Coulomb correction to an arbitrary matrix vertex — ir(p, e; q, w), corresponding to the 
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process, shown in Fig. 25: 



/duj' da' 
^ (2^ ^(q'. ^') G(p + q', e - w') r(p + q', e + a;'; q, lu) G{p + q + q' ,e + uj + oj') 



2tt (27r 



27r (27r)2 gr' g„g' + ^(wg')2 + ,^2 

e-i-cu + v(p + q') -"E ^, , , . e - lu - izu + v(p + q + q') ■ "E 
7 — : ^ TT. TTTT r(p + q , e + w ; q, w) -, — -. -. rr^ ^^ -7 

^maa: CO 



,, , , VFT^ w^T-e^r-E/2 

2nv J "^ J '''^ g^+ VFT^ (^' + ^')' 

4.9 . e™ax ? , (l-sinV)r-sin2(^SrS/2 
m dip ■ — 



7r2A/' £.min J I + g sm Lp 





8 6r 



/(5)r-/(5)(r + srs/2) 



(D2) 



Here we have assumed that the dependence of r(p, e; q, uj) on p and e is weak. If this dependence conies entirely from 
the renormahzation, which is true in our case, its weakness is due to the smahness of the parameter 1/J\f. 
Let us repeat the calculation of the self-energy (50) due to electron-phonon interaction: 



where the phonon mode index /i runs over the two modes belonging to the E2 representation and the two modes 
belonging to the E[ representation, the corresponding matrices (AE)^ being — A^Ey, A^E^,, A^^z, AyJ^^. Note that 
the here the integral is dominated by nj ^ max{e, a;^}, in contrast with the Coulomb self-energy where we had zu ^ vq. 
Thus, to calculate the integral in Eq. (D3) we simply approximated 

D^{q,iw)^ ^I^^-2tt5{w), |tn|>w^. (D4) 

The sum of the two diagrams in Fig. 24 gives the following expression for the effective two-electron vertex [we 
denoted p = (p, e), q = (q, uj) for compactness]: 

ep(2), , , ■ f A' Y- ^' V27a2 _F^V27a2 

^^ (2'2'2) = V(^g 2Af^, 4 2MU.,, 4 D,iq_-q_)DAg)x 

X (AE)^G(£ + g') (AE)^, ® [(AE)^G(£' - g') (AE)^, + (AE)^, G(£' - g + g') (AE)^] , (D5) 
Due to the condition (D4), wc can approximate 

G{p + g') « -G(£' - q') « G(£' - g + g') « -^^ ■ (D6) 

Then the frequency integral is simply 

dio' , ,, , ,, 2(0;,, +0;,,/) , ^ , , 

-D,i.- .') D, (.') ^ ^._;;+c.;'.o)2 - ^-+-' ('^) • (^^) 

Taking the q'-integral, we obtain 

<5r(2)(p,p';g) = !^in|^^ ^ |1 (AE)^E,(AE)^, ® [-(AE)^E,(AE)^, + (AE)^,E,(AE)^]1 . (D8) 

The matrix expression appearing in the braces is evaluated for each of the 16 combinations of the indices /i, /i' in 
Table IV. Adding up the contributions, we obtain Eq. (134b). 
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H,fi' 


-A.Sj, 


A^E^ 


A:.E. 


AyE, 


-A,S^ 








AyE^ (g) AyE^ 


A,E, ®A,E, 


A,S, 








AyE^ g)AyE^ 


Aa;E^ ® A^E^ 


A.E. 


AyEz ® AySz 


AyT.z ® AyEz 





A^Ey (g) AzEy + A^Eo, (g A^E^z, 


Aj/S^ 


Aj^E^ (g) Aj:Ez 


A^^E^ (g) AajEz 


A^Ey (g) A^Ey + A^E^ (g A^E:^, 






TABLE IV: The matrix expression in the braces in Eq. (D8), evaluated for all 16 combinations of the phonon indices /J., n' ■ 



APPENDIX E: POLARIZATION OPERATOR 

Polarization operator with arbitrary matrix vertices S^ , S^ , i, j = x,y, z,0 (we denoted Eq = 1) i can be conveniently 
calculated in the coordinate representation:^^ 



Af 



-zHy (r - r', i - t') = — Tr^{S,G(r - r', t - t') SjG(r' -r,t'- t)}, 



{idt + iv-E ■ V) G(r, i) = 1 (5(r) (5(t). 



(El) 
(E2) 



If A matrices are also present in the vertices, tracing of them is trivial: the matrices in the two vertices must coincide 
for the trace not to vanish, then their product is equal to the unit matrix. 

It is convenient to switch to the imaginary time r — it, idt ^ —dr- Then the Green's function in the coordinate 
representation can be found by using the analogy with the 3D Coulomb problem. Namely, we introduce the third 
dimension z = vt, and notice that 



i-dr + ivT, ■ W){-dr - ivi: ■ V) = v^WJo. 



(E3) 



Since the inverse of the 3D Laplacian is the Coulomb potential, we obtain 

1 WT + zS • r 



G(r, r) = {dr + wS • V) 



Attv^/v^t^ + r2 4^(w2r2 + r^f/^ 



(E4) 



Using auxiliary relations (here i,j — x,y, z), 



Tr^lSiSfeSjS;} = 4:SikSji + AS^iS^k - 4%4( 



did J 



1 SxiXj — 25ijR^ 



i?2 i?6 

we calculate the polarization operator (here i,j = x,y): 



, R^ =x^ +y 






(E5a) 
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(E6a) 
(E6b) 
(E6c) 
(E6d) 



Using the 3D Fourier transform 



/ 



d-'R 



oiQR. 



(fn n 



327r2i?2 J (27r)3 4i?2 



16 g' 



(E7) 
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we obtain n(q, lu) (up to a q, a;-independent constant coming from r,T ^> 0): 



noo(q,w) == - 

no*(q,^) = - 

ny(q,tj) == - 



? 9 



16i;2 ^i;2q2 _ ^2 ' 

gi^j A/" w^ 



^2 16i;2 ^^;2^2 _ ^2 

2A/' 



(*.-t)i^^^^v^. 



n2^(q,tj) = — -^ Vw2g2_t^2^ 

16?''^ 
n,,(q,L^) - 0. 



(E8a) 
(E8b) 

(E8c) 

(E8d) 

(E8e) 
(E8f) 



Given the polarization operator, one can find corrections to the optical phonon frequencies due the electron-phonon 
interaction: 



A 



K „,2 



K 



D],\q,u;)^^v'U,_,ici,Lu) = => 5u: K{ci) ^ ^ \ v\^ ^ u^l, 



2 

^ „,2 

2 



4 
Ar 



D^\{(1,u)-^v^IIt{^,oj) = Q => 5a>r,L(q)«^Vt;2^2_^2 



F' 



D^}^{q,Lu) -^v^ULiq,Lu) = ^ (5cJr,T(q) ~ -^^ 



8 yjv^q^ — ojp 



(E9a) 
(E9b) 
(E9c) 



At wg < ^JJ^J_ the square roots are imaginary which corresponds to the decay of phonons into the continuum of electron- 
hole pairs. 
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